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Abstract 

We compute the amplitudes for the two-body decay of B mesons into longitudinally 
polarized light vector mesons at next-to-leading order in QCD. We give the explicit 
expressions in QCD factorization for all 34 transitions of a heavy-light B meson into 
a pair of longitudinal vector mesons p, lo, 4>, K* within the Standard Model. Decay 
rates and CP asymmetries are discussed in detail and compared with available data. 
Exploiting the fact that QCD penguins are systematically smaller for vector mesons 
in comparison to pseudoscalars in the final state, we investigate several methods 
to achieve high-precision determinations of CKM parameters and New Physics 
tests. We propose a method to use V-spin symmetry and data on B^ K*j^K*^ 
to constrain the penguin contribution in B^ P^Pl- CP violation in the latter 
decay together with a measurement of sin 2/3 determines the unitarity triangle with 
high accuracy. We show that CP violation in B^ — > p^p^ and B^i ipKs alone 
implies \Vub\ = (3.54 it 0.17) • 10~^, presently the most accurate determination of 
this quantity. 



PACS: 12.15.Hh; 12.39.St; 13.25.Hw 



1 Introduction 



Charmless hadronic S-meson decays provide us with important information on the 
physics of flavour, within and beyond the Standard Model. After the very successful 
first generation of 5-factory experiments, BaBar and Belle, the interest in this field is 
re-inforced by the upcoming start of the LHC with the dedicated LHCb experiment [1]. 
In the longer run the second generation S-factory projects SuperB [2] and SuperKEKB 
[3] promise excellent future opportunities [4,5], as well as the LHCb upgrade options [1]. 

In the present paper we investigate in detail the class of hadronic two-body decays 
of B^, Bd and Bg mesons with longitudinally polarized light vector mesons (p^, a;, 
0, K*^, K*^, K*^) in the final state. These can be computed systematically using QCD 
factorization in the heavy-quark limit [6,7,8]. In comparison, B decays into vector mesons 
with transverse polarization are suppressed by powers of KqcD/mb- Their amplitudes 
do not factorize, as indicated by infrared singularities in convolution integrals, and this 
introduces a model-dependence in the computation of the related obscrvables [9,10,11]. 
For this reason we do not include final states with transversely polarized vector mesons in 
our present analysis, although they raise issues of interest in their own right [9] . Rather, 
we will pursue the physics of B ^ VlVl decays, which comprise 28 calculable channels 
and 6 annihilation modes within the Standard Model and offer a rich phenomenology 
by themselves. The various decay modes probe different types of amplitudes, such as 
tree, QCD and electroweak penguin, or weak annihilation. The large number of channels 
will allow us to test our understanding of QCD in hadronic B decays and at the same 
time to extract fiavour parameters and search for New Physics effects. An important 
motivation for the study of final states with two vector mesons is the considerably smaller 
size of QCD penguin contributions, as compared to the decay into pseudoscalars. This 
improves our control of QCD effects in several observables of interest for fiavour physics. 
Mixing- induced CP violation in B^ PlPl an important example. 

Several studies of S — > VV decays based on factorization in the heavy-quark limit 
have been made in the literature [9-18] . So far the only comprehensive analysis, including 
all channels, has been given in [11], where also transverse polarization is discussed. The 
other articles have addressed various aspects of i? — > VV decays, investigating particular 
channels, polarization effects and, in part. New Physics contributions. 

In this paper we present a detailed analysis of the complete set of 5 — > VlVl decays 
within the Standard Model. We work to next-to-leading order (NLO) in QCD, using 
factorization in the heavy-quark limit. Power corrections to this limit, in particular from 
weak annihilation, are estimated using a model description [6]. Our presentation includes 
some additional aspects and several new applications of S — > VlVl transitions. We give 
detailed analytical expressions for all 34 S — > VlVl decay amphtudes at next-to-leading 
order. We suggest an alternative treatment of long-distance electromagnetic penguins, 
which contribute in channels with the neutral vector mesons p°, u; or 0. The effects of uj-(j) 
mixing are briefly discussed. We apply the results in new ways to the phenomenology of 
the unitarity triangle: we use the formulation of [19,20] for the analysis of CP violation 
in Bd — > PlPI-i we propose an approach to independently constrain the B^ — > PlPl~ 
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penguin with K'^K^ and V-spin symmetry, and we extract accurate values of 

the CKM quantities p, fj, 7, a and \ Vub\- 

The paper is organized as follows. Section 2 collects basic ingredients of the calcula- 
tion, in particular the effective Hamiltonians, form factors, light-cone wave functions and 
meson projectors. Section 3 presents the NLO results for the B — > VlVl decay ampli- 
tudes. Input parameters and experimental results are summarized in section 4. Section 
5 contains the phenomenological analysis, with a discussion of branching fractions and 
CP asymmetries and with precision determinations of the unitarity triangle. A brief 
comparison with the literature is given in section 6, before we conclude in section 7. Our 
treatment of long-distance electromagnetic penguins is discussed in the appendix. 

2 Preliminaries 

The effective weak Hamiltonian for charmless hadronic B decays, without change in 
strangeness {AS = 0), is given by [21] 

7if/=° = ^ ^AJcig? + C2gf+ QQi + C^^Q^^ + CsgQsg) +h.c. (1) 

^ p=u,c ^ i=3,...,10 ^ 

where the elements of the CKM quark-mixing matrix V enter as Xp — VpbV*^, Ci are 
Wilson coefficients, and the operators Qi read 





= {pb)v- 


-A{dp)v-A , 


0^2 


= (Pibj)v- 


-A{djPi)v-A , 


Qs 


= {db)v- 


-A^q {qq)v-A, 


Qa 


= {dibj)v- 


-A {Qj(li)v-A , 


Q5 


= {db)v- 


-A {m)v+A , 


Qe 


= (dibj)v- 


-A ^q iqjqi)v+A , 


Qv 


= {db)v- 


-A , 




= {dibj)v- 




Q9 


= {db)v- 


-AYqleq(qq)v-A, 


Qw 


= {dibj)v- 


-A \^q{qjqi)v-A 








9 

^ 





(2) 



Here i,j are colour indices, are the quark charges, and the sums extend over q — 
u,d,s,c,b. The Wilson coefficients Cj will be taken at next-to-leading order (NLO), 

using the treatment of electroweak contributions described in detail in [6]. The sign 
conventions for the electromagnetic and strong coupling correspond to the covariant 
derivative = dn + ieQfA^ + igT"'A^^. With these definitions the coefficients Cy^, 
Csg are negative in the Standard Model, which is the convention usually adopted in the 
literature. 

The effective Hamiltonian for charmless decays of B mesons with AS — 1 can be 
obtained from (1) by interchanging d- and s-quark labels. The CKM factors governing 
these transitions are then = Vpiy*g. 
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To obtain the amplitudes for B — > VilV2l decays from the Hamiltonian, the matrix 
elements of the operators Qi have to be computed in QCD factorization [6,7,8]. To 
lowest order the matrix elements are expressed in terms oi B ^ V form factors and 
vector- meson decay constants. The required form factors are defined by (see e.g. [22]) 



{V{p,e)\q^^^^h\B{pB)) = 2mvAo{q^f-^q^+{mB + mv)A,{q') 



-M{,q 

VilB + my 



\pB+Pr-'^^^q^ 



(3) 



e)\qrh\B{pB)) = -2i e^^'"^pBuPpe„ (4) 

mB + Triy 

where the momentum transfer is q — pb — p and the totally antisymmetric tensor ei^'^P"' 
is normalized by e^^^^ = — 1. 

The vector-meson decay constant fy is given by 

{V{q,rj)\q-ff,q'\0) ^ -ifymyrji^ (5) 

for a vector meson with flavour content V = q'q. The corresponding matrix element 
where 7^ 7^75 is zero. For energetic vector mesons V{p, e) with longitudinal polariza- 
tion 

e'-— (6) 
my 

up to corrections of second order in my/niB- 

The factorized matrix element of a (l^ — ^4) ® — A) operator then reads to lowest 
order (0;°) 

{V,LV2L\{qib)y-A{q2q2)v-A\Bg) = imlA^^'''{m\)fy, (7) 

where Vi = qqi and V2 = ^2^2 (^i 7^ 92)- 

The corrections at higher order in are expressed in terms of calculable hard- 
scattering kernels and meson light-cone distribution amplitudes. The latter quantities 
enter through the meson projectors in momentum space. For the B meson this projector 
is given by 

bq = ^(^B + ruBh, [0si(O + 0S2(O fi] (8) 

In our notation (bq) denotes a matrix in Dirac space displaying the flavour composition 
of a Bq meson in the initial state. The 4-vector = (1,0,0,-1) is chosen to be in 
the direction of the recoiling meson Vi. The parameter ^ is the light-cone momentum 
fraction of the spectator quark q. The distribution amplitudes are normalized as 

-1 rl 



I 



f dCMO-h [ dC<f>B2{0-0 (9) 

Jo 



In the present analysis does not enter the results and (pBi appears only through the 
first inverse moment 

d^^^^ (10) 



f 

Jo 
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which defines the hadronic parameter = 0{Aqcd)- Colour indices have been sup- 
pressed in writing (8). They are taken into account by replacing bq biQj and including 
a factor of Sij/Nc on the right-hand side. 

For a longitudinally polarized vector meson in the final state with flavour content 
^2?i and momentum p, the projector can be written as [23] 

^2?! = ^ - — my^— (11) 

Here x is the momentum fraction of the flnal-state quark qi and 

P 

E = xp^ + E + -^n" (12) 
2xmB 

P 

kit = xpf" -E + ^^n^" (13) 

(with X — 1 — x) are the momenta of Qi and q2, respectively, is a light-like vector 
with spatial direction opposite to p^: If p ~ n+ (p ~ n_) then n — {n — n+), where 
< = (1,0,0,±1). 

The function 0||(.x) is the light-cone distribution amplitude of leading twist for a 
longitudinal vector meson. The subleading-twist amplitude ^v{x) has been treated in 
(11) in the Wandzura-Wilczek approximation. It gives rise to contributions suppressed 
by one power of tvQCDj'^B- We will nevertheless include it in order to estimate the 
impact of this particular source of power corrections to factorization in the heavy-quark 
limit. 

The functions 0|| and $^ can be expanded in terms of Gegenbauer and Legendre 
polynomials, respectively, 

oo 

(\>\\{x) = ^xx ^ a„C^/2(2x - 1) (14) 

n=0 

oo 

^,{x) = 3^a„^P„+i(2a:- 1) (15) 



n=0 



where ao = o^o± = 1- In the Wandzura-Wilczek approximation $^ can be expressed in 
terms of the twist- 2 wave function of a transversely polarized vector meson, (f)±, as 

X = / du^-^- / du^-^^ (16 
Jo u u 



0_L has an expansion similar to (14) and this leads to (15). 
Note that (16) implies 

dx^y{x)^0 (17) 

'0 



/' 

^0 
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even though is not necessarily antisymmetric under x ^ x for general an±- The 

normalization of (p\\^± is dx 4>\\^±{x) ~ 1. 

For phenomenological applications we shall truncate the expansions of 0|| and and 
use, for a particular meson V, 

0f[(a;) = 6xx [l + 3{2x - 1) + 6{5x'^ - 5x + 1)] (18) 

C(^) = 3(2x - 1) (19) 

Taking the vacuum-to-meson matrix element of a local current, the projector (11) 
reproduces (5), (6) 



(^(p)|?i7m?2|0) = - / dxtrj^q2qi = -ifvPn 
Jo 



(20) 



and 

{V{p)\q,Tq2\0)^0 (21) 

for r = 1, 75, 7^75. 

3 QCD factorization in B VlVl decays 

The amphtudes for the AS — decay of a S meson into a pair of light vector mesons with 
longitudinal polarization can be conveniently expressed as follows (the case of AS — 1 
is obtained by replacing d s): 

{v^LV2L\n^^'='\B) -^J2^p {VilV2l\t,' + t;-''^|b) (22) 



p=u,c 



where 



Tp^ = ai(ViV2) 5pu (ub)v-A ® (du)v~A 
+ 02(141^2) Spu {dh)v-A ® {uu)v-A 

+ 03(^1^) {db)v-A ® {qq)v-A 

+ al{ViV2) J2i (qb)v-A ® {dq)v-A 

+ 05(^1 1^2) {db)v-A {qq)v+A 

+ 07(^1^2) idb)v-A ® ^eg{qq)v+A 
+ 09(14^2) {db)v-A leg{qq)v-A 

+ <o(^i^2) iqb)v-A ® leg{dq)v-A (23) 
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Here the summation is over q = u,d, s. The symbol ® indicates that the matrix elements 
of the operators in are to be evaluated in factorized form [6]. The factorization 
coefficients include hard QCD corrections to the S-decay matrix elements at NLO, as 
well as electroweak effects in the systematic approximation of [6]. Note that structures 
with scalar and pseudoscalar currents are absent in (23), in contrast to the case of 
B — > Kt: considered in [6]. Because (pseudo)scalar currents cannot create a vector 
meson from the vacuum, these structures can give no contribution to 5 — > VV decays. 

The term T^^nn.d ^22) describes the effects of weak annihilation. These are power 
suppressed in the heavy-quark limit and cannot be computed in QCD factorization. We 
shall use model calculations to estimate this important class of power corrections to the 
leading, factorizable amplitudes. Weak annihilation will be discussed in Section 3.3. 

3.1 Results for the parameters ai 

The factorization coefficients can be written as = + ajji. We find 



= Ci + 



02,1 



C2 + 



03,1 = C's + 



= C4 + 



C2 
Nc 

Nc 



1 + 



1 + 
1 + 
1 + 



Cpaa 
An 

CpCXs 

An 
An 



V 



V 



V 



C5 CpOts V-irl. 

ryVy 



Nc An 



C2 CpT^as rj 
C4 Cpnas rr 



V1V2 



«4,n 



C3 CpTias 



V1V2 



'10,1 



_ ^ C'lo 



1 + 
1 + 



An 

CpOLs 

An 



1 + 

An 



, " pp,EW 



9n 



a 



An 



^9,11 



Nc An 



V 



9nN, 



Nc Nc 

Cg CpT^OLs 



Nc Nc 

Cio CpT^Ois 



OlO,II 



Cg CpHQ.s 



Nc Nc 



V1V2 



V1V2 



(24) 



where Ci = Ci{fi), as = Q;s(/i), Cp = {N^ — l)/{2Nc), and Nc = 3. The hadronic 



quantities Vy \ Hy^y^, Py2, Pv,zi PvS"" ^v^"^ 1 ^"^^ ^vln"" given below. All indices 



jP,EW pp,EW 



,p,EW 
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Figure 1: Vertex diagrams. 



V in Vv, -Py, are understood to refer to the emitted meson V2. 

Contributions suppressed by one power of ^-Qco/'^nb that arise from the twist-3 com- 
ponent of the vector-meson wave function have been included in the above expressions. 
They are related to the scalar penguin operator {qb)s-p{sq)s+p and come with a factor 



y. , 2mvf^ifi) 2my/^(lGeV) 



mb{n)fv 



mb{mb)fv 



asinib) 



a,(lGeV) 



Cf/Po 



(25) 



Here mb{fi) is the M5'-mass of the b quark at scale fi, and /3o = 23/3 for / = 5 flavours 
of quarks. 



3.1.1 Vertex and penguin contributions 

The vertex corrections (Fig. 1) are given by 



12 In 


nib 




1^ 


12 In 


rUb 




fi 


„ /I 


- 2x 



1 — X 



■\nx — in 



+ 



2L2(x) -ln^x + 



2 Inx 
1 — X 



— (3 + 2z7r) In a; — (x 1 — x) 



= [ dx [2L2(x)-ln2x-(l + 2i7r)lnx-(x^l-x)] $^(x) 
Jo 

where L2{x) is the dilogarithm 



Lo(x) = - / dt 



t 



The expansion of 0j| in Gegenbauer polynomials gives 

dxg{x) (p\{x) = -\- 3«7r + - ?>iTX^ "1^ - ^ «^ + • • • 



(26) 
(27) 

(28) 
(29) 




Figure 2: Penguin diagrams. 



Replacing g{x) hj g{l — x) leads to a change of sign in front of the odd Gegenbauer 
coefficients on the right-hand side. 

Next, the penguin contributions (Fig. 2) are 



pp 



4: m, 2 

-In — + - - Gv{sp) 

6 jJ, 6 



^ln^ + ^-Gv(0)-GHl) 



+ (C4 + Ce) 

"-^ dx 



^ln!I^-3Gy(0)-Gv.(..)-Gy(l) 

6 fJ, 



-2C: 



off 



lo 1 - a; 
pP'^^w ^ ^ ^^^^ 



X] 



4, rn.2 



dx 
1 — X 



(30) 



where = and Sc = {mc/mby. Small contributions from the electroweak coefficients 
C7, . . . , Cio are consistently neglected in within our approximation scheme. Also, 
the very small corrections from C3, . . . ,Cq in (and in Py^^ , see (35) below) are 

omitted for simplicity. 
The function G'y(s) is 



Gvis) = I dx G{s — ie,l — x) (f)^^ (x) 



G{s,x) = —A / duu{l — u)\n[s — u{l — u)x] 



(31) 



2(12s + 5x - 3a;lns) A^/As - x {2s + x) 



9x 



3x3/2 



arctan 



X 



4:S — X 



(32) 
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Expanding in Gegenbauer moments one finds 



Gvisc 



V 



a 



V 



^-|lns.+ ^ + ^ + 3(8 + 9ar + 9ar). 



c 



+ 2(8 + 63q;[ + 2Ual)sl - 24(9q;]; + mal)sl + 2880q;2 s; 
1 + 2sc + 6(4 + 27a];' + ISa^sl 



- 2>&{9a\ + l^aX)sl + 4320a2^s 



.y„4 



(2 arctanh\/l — 4sc — 



+ 12s^ 



4 



1 + 3al' + Qal - - (l + 9a]f + 36a^) 



+ 18(a[ + 10a^)s^ - 240^2'' s. 



,y„3 



(2 arctanhVl — 4sc — ^tt) + . . . , 



2«7r 



V 
"2 



G'y(O) 

G'y(l) = ^ - 6^3 TT + ^ - ( ^ - 36^3 TT + 12r 



85 



47r2 



155 



+ _ 504 V3 TT + 1367r^ ] + . . . 



(33) 



The function Gv'('S) and its expansion in Gegenbauer moments liave tlic same form as 
Gk{s) in tlie case of 5 — > Kn discussed in [6]. Likewise tlie integrals proportional to 
Cy^ and C|J in (30) are similar to those in S — > Xtt. They read 



dx 



lV , 



x) = 3(1 + < + + . . . ) 



X 



(34) 



The twist-3 terms from the penguin diagrams are obtained from the twist-2 terms 



by the replacement 



X) 



~(cff 



$^ (,x) , except for the terms proportional to Cf^ and C^g . 



Here the factor of (1 — a;) in the denominator of the integral in (30) is canceled by the 
twist-3 projection. An important difference between the twist-3 penguin contributions in 
B ^VlVl and B — > Xtt arises from the different properties of the twist-3 wave functions 
in these two cases. Since Jq dx^^{x) = it follows that the contributions from C^^ and 
CgJ vanish in the former case. The same holds for all x-indepcndcnt constants in the 
hard-scattering kernel, in particular for the scale and scheme dependent terms. We then 
find 



DP _ 

-^ys 



r, Gv{sp) + C3 (G'v(O) + Gv{l)) + (C4 + Ce) (3(5^(0) + Gv{sc) + G'v(l)) 

-(Ci + iVeC2)G'v(sp) (35) 



pP,EW 

-^ya 
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with 



(36) 



Gv{s)= / dxG{s-ie,l-x)^^{x) 
Jo 

Using the asymptotic form of the wave function ^y{x) — 3{2x — 1) leads to 
Gv{sc) = 1 — 36sc + 12sc\/l — (2 arctanh\/l — 4sc — "^tt) 
— 12s^ (2 arctanhVl — — ot)^ , 

G'y(0) = l, Gv(l) = ^7r2 + 4y37r-35 (37) 



Finally, we give the electromagnetic penguin contributions {p = u, c). For 

intermediate charm, p — c, these are calculable in perturbation theory and read 



A , nih 2 4 , mc 

-In— H h-ln — 

3 /X 3 3 TTib 



(38) 



In the case of the up-quark loop, p = u, the amplitude becomes sensitive to additional 
long-distance dynamics, which is not strictly calculable. Using a suitable hadronic rep- 
resentation of the light-quark loop, we estimate 



3U,EW 
V,n 



(Ci + N,C2) 



rub 10 47r^ v-^ 
- m 1 > 

9, II Q S ^ 



n 



2% m?/ 2 , m 
-I + -ln — 



3 tr 



+ 



9 

2 t. 



m 



r=p,oj 



V 



m 



V 



m 



V 



(39) 



where tc = 47r^(/^ -|- f^). This point is discussed further in appendix A. 

The authors of [23] factorizc the term (39) into a short-distance and a long-distance 
part, separated by a scale u. The short- distance part is equivalent to (38) with rric 
replaced by u. The long-distance part is not considered explicitly in [23]. Our treatment 
is consistent with the framework of [23], but supplies a concrete model representation 
for the long-distance contribution of the electromagnetic penguin. 



3.1.2 Hcird spectator scattering 

The hard spectator interactions (Fig. 3) determining the coefficients Uiji in (24) are 
governed by the quantities 



HviV2 — 



^2 /iS-^Vi 



(0) 



(x) J (y) + rl^ (y) (40) 



TT' _ 



Vi 



^(0) 



^ dx 



dy 



<l>ii'iy) + rl'^^v'iy) (41) 
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Figure 3: Hard spectator diagrams. 



Strictly speaking, only the twist-2 components of the vector meson distribution ampli- 
tudes contribute at leading power and are consistently calculable in the present approach. 
The twist-3 part described by leads to terms with logarithmic endpoint singularities, 
but these terms are suppressed by one power of A/mfe. We shall include them in our anal- 
ysis using a simple model, in order to estimate the potential impact of power corrections 
from this source. Following [6], we parametrize the endpoint singularity by 



The logarithm comes from cutting off the lower range of integration at ymin = ^h/i^B, 
and ph, 4>h are real model parameters to allow for a complex Xh and a deviation from 
the default value \n{mB/Ah). The integral over in (40), (41) then becomes 



Throughout we use ph = V-^h with Ah = 0.5 GeV as the scale in the spectator- 
scattering contributions. 

3.2 B — ^ V\Ly2L decay amplitudes 

The transition operators T*^, describe a total of 28 two-body decays of , and 
Bg into the charmless vector mesons K*^ , K*~ , K*^, K*^, p"*", p~, p°, u and 0. There 
are 15 AS* = {b ^ d) and 13 AS" = 1 {b ^ s) transitions. In this section we give the 
expressions for the amplitudes of these processes in terms of the factorization coefficients 
Oj. All light vector mesons are taken to be longitudinally polarized, that is, p~p° here 
means p2Pl- We use the abreviation 




(42) 




(43) 



Av,v2 = i^mB^o 




(44) 
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suppressing the dependence of Ay^vb on the i?-meson flavour in the notation. 

The decay amplitudes, up to the factor (44), are conveniently obtained from the 
transition operator by the following 'bosonization' of the bilinear quark currents in (23): 



{ub)v-A 

{db)v-A 

{sb)v-A 
{du)v = p~ 

{ds)y = K*\ 



+ + BsK* 



v2 

B-K*- + BdK*'' + Bscf) 

{ud)v = {sd)v = 

{su)v = K*-, {us)v = K*+ 



[uu)v = 



^/2 ' 



{dd)v = 



V2 ' 



{ss)v = 



(45) 

(46) 
(47) 

(48) 



In our notation the charge of B mesons (hght mesons) corresponds to that of particles 
in the initial (final) state. Note that the axial vector parts of the light-quark currents 
don't contribute for final-state vector mesons. Insertion of these expressions in T'^ (T*) 
from eq. (23) generates all AS* = (AS* = 1) amplitudes. For a specific process 
B — > V1V2 the amphtude is found as the coefficient of (5^1)^2 (and of {BV2)Vi if 
Vi 7^ ^^2). This procedure automatically keeps track of all sign and Clebsch-Gordan 
factors. Note, however, that an extra symmetry factor of 2 has to be included for 
amplitudes with two identical particles in the final state. From the structure of T'^^^^ it 
follows that the coefficients 03, and a^, Og always appear in the combination 03 + 05 
and + ag, respectively. Representative numerical values for the coefficients can be 
found in appendix B. 

The AjS = transition amplitudes then read (a summation over p — c is under- 
stood): 



\/2A{B- 



rp') 



Xu{ai + 02) + -Ap(af + al + afo) 



(49) 



V2A{B- p-uj) = [Kai + \p{al + a?o)] A 
+ 



A„a2 + \p {^l + 2(03 05) + ^(a? + «9 - «io)^ 
A{B- p-(t>) = Ap 
A{B~ K*-K*^) = Xp 



1 



03 + 05 - 2(^7 + ^9) 



"4 2 10 



Ak*k* 



A{Ba ^ pV°) = 



A„a2 + Ap I - -(a? a^) - -afo 



(50) 

(51) 
(52) 

(53) 
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Ka2 + Ap -al + -(a? + ag) + -a' 



10 



Aua2 + Xp ^af + 2(a3 + 05) + ^(o? + Og - Oio)^ 



pU! 



A„a2 + Ap ( of + 2(03 + as) + ^C^? + «9 " «io) 



«3 + 05 - ^(a? + al) 



V2A{Ba ^ pV) = -Ap 

V2A{Bd cu(t>) = Ap 

A{Bd p^p~) = [A„ai + Ap(a^ + afo)] >lpp 
A{Bd = A 



03 + 05-^(07 + 09) A 



A„a2 + Ap ( -al + ^(o? + af) + ^afo 



V2^(5^ ^ /sT -tu 
^(^3 ^ i^*V) = Ap 



Aua2 + Ap ( + 2(03 + a^) + ^("t + «9 ~ «io) 



p p 
04 — 2'^io 



^fliftT* + Ap 

A(Bs K*+p-) = [A„ai + Ap(a^ + a?o)] ^/^v 



03 + 05 - -(o? + a^) 



A 



K*p 
K*4> 



The amplitudes for — 1 transitions are found to be: 



^A{B- ^ i^*-pO) = [A;ai + A;« + a?o)] A^k* + 



A>2 + -A;(a? + al) 



V2AiB- ^ K*-oj) = [A>i + A;« + <o)] 



+ 



Al 



a>2 + a; 2(03 + 05) + -K + 



04 + 03 + 05 - -(a? + + afo) 



A{B- ^ x>-) = a: 



04 — 2^10 



ApK* 



13 



-a 



10 



(68) 



A>2 + A;(^2(a3 + a5) + i(a? + a^)^ 



+A' 



iuiK* 



04 + «3 + 05 - 2 (^7 + + ^lo) 



A{Ba ^ = a; 

A{B, ^ = [X'^a, + A;« + afo)] A,k 



A>2 + :^A;(a? + ) 



A 



4>p 



A{B, 



X'^a2 + Ap I 2(a3 + a^) + - (a? + af) 



A, 



2A' 



04 + 03 + «5 - :^(«7 + 09 + afo) 



A, 



A(B, ^ = [A>i + X'^ial + <o)] 



a 



10 



Ak*k* 



(69) 

(70) 

(71) 

(72) 

(73) 

(74) 
(75) 
(76) 



3.3 Weak annihilation amplitudes 



The decay mechanism of weak annihilation (Fig. 4) gives contributions to the amphtudes 
for B VlVl decays that are supressed by AQcn/nT-b- These power corrections are 
not calculable in the usual factorization framework. This is indicated by end-point 
singularities from the integrals over light-cone momentum fractions in a hard-scattering 
ansatz. We shall use the model of [6], which is based on this ansatz together with a 
cut-off procedure, to estimate the impact of annihilation effects on the leading decay 
amplitudes. Following the notation of [6] we write 



\B) 



(77) 



p=u,c 



with Xp = VpbV*^ and 



T 

p 



ann,d 



= Sup {Srd bl + 5ru b2 (Trf ) 63 0"^ + Srd &4 tr((j) 

+ ^^re,a^ + ^5,,6rtr(ga) 



(78) 
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Figure 4: Annihilation diagrams. 



Here the index r = u, d, s denotes the flavour of the spectator quark in the B meson 
and Q = diag(2/3, —1/3, —1/3). The corresponding formulas for A^* = 1 transitions are 
obtained by interchanging the labels ci ^ s in the expressions (77), (78) for A^* = 0. In 
particular, Xp is then replaced by = VpbV*^. 
The operators 

<= E (?"92)x(gV) (79) 

q'=u,d,s 

encode the valence quarks of the final state mesons. Matrix elements of the product of 
currents in (79) are defined as 

{V,lV2l\ji X J2\B) = icfBfvJv, (80) 

where c = 0, 1, ±l/\/2, etc., is the appropriate Clebsch-Gordan coefficient and symmetry 
factor, relating the currents ji and j2 to the mesons Vi, V2. The coefficients bi are given 
by 



c 



C;A'^ + C,{^^ + Ai) + N,C,Ai 



LEW 



C,A', + C^{A'^ + Ai) + N,C,Ai 



(81) 
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and correspond to current-current annihilation (61, 62), penguin annihilation (63, 64), and 
electroweak penguin annihilation (b^^ , 6f These coefficients depend on the final-state 
mesons, bi — 61(^1 V2), but this dependence will be left implicit in the following. Finally, 
the A^^^ read 





— na 


A{ 


= 0, 




= Tia 


A{ 


= 0, 


^3 


= ttq; 


4 


= ttq; 



^0 



+ 



1/(1 - yx) yx^ 



yx ) 



dxdy <; 0[[i(|/)0|*|'2(a; 



x{l — yx) y'^x 



V2, 



yx 



^nas / dxdy \ (l)X^{y)^l^{x)r 



.V2 



2x 



yx{l — yx) II 



2^ 



yx{l - yx)}' 



^0 I y 



2(1 + a:) 
yxP' 



(82) 



The superscript i (/) denotes gluon emission from the initial- (final-)state quarks, as 
shown in Fig. 4 (c) and (d) ((a) and (b)). The subscript k indicates the Dirac structure 
of the four-quark operators, Vi = {V - A) ® {V - A) {k = I), {V - A) ® {V + A) 
{k^2), {-2){S-P) ®{S + P) {k = 3). The various quantities in (81) wiU be evaluated 
at the scale /x^ = similarly to the spectator-interaction terms. 

For the numerical estimate of weak annihilation the IR-divergent quantities A^j^-^ in 
(82) will be parametrized by 



Xa 



TUB 



— = (l + pAe'^^)ln , 
X ^ 'Ah 



(83) 



The quantity Xa is the cut-off regulated integral Jj^^^^^ dx/x = \n(mB/ A^), with scale 
Afi = 0.5 GeV, modified by a phenomenological magnitude pa and phase (j)A [6]. Using 
SU{3) flavour symmetry and the asymptotic forms of the meson wavefunctions 0|| and 
one finds A\ = , = and 



lS{XA-^ + -]+18{rl)'{XA-2)' 



A\ Ki na 
Ai fti -SGnasrl {2X1 - + 2) 



(84) 



We next give the results for the annihilation amplitudes in terms of the coefficients 
hi, where we define 

BviV2 = fBfvifv2 (85) 
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The decay constant Jb depends on the flavour of the decaying B meson, even though 
this is not made exphcit in the notation for the -ByiVj- The following expressions can 
be efficiently obtained with a procedure similar to the one described at the beginning of 
sec. 3.2. Typical numerical values for the coefficients bi are given in appendix B. For 
the channels with A^* — the annihilation contributions read 



A„n(5- ^ p-p°) = 

\/2Ann(5" ^ p-u;) = [\u 2b2 + (A„ + A,) 2(63 + 63"^)] 5 

Ann(5~ ^ p-(t)) = 

A,UB- ^ K*-K*^) = [A„ 62 + (A„ + Ae) (63 + Bk*k 



Xu bi + {Xu + Ae) (^63 + 264 - + ^^r) 

A.6i + (A„ + Ae) (-63 + ifefw + 3^EW^ 
A„ b^ + (A„ + Ae) (^63 + 264 - + ^^r) 



5 



pa) 



A^n{Bd ^ p+p-) = A„ 61 + (A„ + Ae) (^bs + 264 - ^fef"^ + ^^r^ 

Ann(5, ^ = (A„ + Ae) (^63 + 2^4 - " b^^^ B^^K* 

V2Ann(S, ^ K*V') = (A„ + Ac) (^-63 + ^i^V 
V^^a.n(S, ^ K*''^) = (A„ + Ae) (^63 - ^&r) ^i^*.- 
Ann(^, ^ ir*V) = (A„ + Ae) (^63 - ^^^^ ^i^*^ 
Ann(S, ^ X*+p-) = (A„ + Ae) (^63 " ^^^^^ Bk 



(86) 
(87) 



(89) 

(90) 

(91) 

(92) 

(93) 
(94) 

(95) 
(96) 

(97) 
(98) 
(99) 
(100) 



In addition there are two AS" = decay modes that proceed only through annihilation 
diagrams: 

Ann(Sd ^ K*+K*-) = 



A„6i + (A„ + Ac) (^264 + ^6r 



(101) 
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>^ann(-Sc; 



(A„ + Ae) (264 - 



(102) 



For the annihilation amphtudes with = 1 we obtain: 

V2 Ann(S- ^ X*-p°) = [a; 62 + (A; + a;) (63 + 6D] Bk*p 
^ i^*-a;) = [Kh2 + (A^ + A^) (63 + &D] 

A„n(i?- - i^*-0) = [a;, 62 + (a; + K) {h + 5/^*0 

^ann(5- ^ i^^V") = [A^, 62 + (A^ + A^) (63 + &D] 
V^Ann(S, ^ X*V°) = + A'J (^-63 + ^6^) 

V2 Ann(S, ^ K*'u;) = (A; + A'J (^h - 

Ann(S, ^ X*V) = (AL + K) (h - ]^T\ Bk*^ 



AUBd - K*-p+) = (A; + K) { h - -hT ] Bk* 



Aa,nn{Bs ^ U;0) = 

= (a:, + A'J (263 + 2h - - &r) b^^ 



"^ann(-^s 
A,nn(Bs ^ X*+X*-) 



5 



(103) 

(104) 
(105) 
(106) 

(107) 

(108) 

(109) 

(110) 

(111) 
(112) 

(113) 

(114) 



Ann(S, ^ ^*°X*0) = (A:, + AJ (^63 + 264 - - (115) 

In the case of AS" 1 transitions there are four pure annihilation modes. Their ampli- 
tudes have the form: 



^ann(-Bs ^ P+P ) = 



a;6i + (a; + aj (264 + -64 



EW 



5, 



PP 



EW 



5, 



PP 



1 



,EW 



B,. 



A;fei + (A; + Aj 



B 



(116) 

(117) 
(118) 
(119) 
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Table 1: Experimental results [24] for CP-averagcd branching ratios and lon- 
gitudinal polarization fractions fi of B ^ VV decays. Here the B meson 
is either a B~ or a B^. The branching ratios for B — > VlVl have been 
obtained as fLB{B VV). Also shown are the direct CP asymmetries 
Acp = -C= {B{B ^ /) - B{B ^ f))/{B{B ^ /) + B{B ^ /)). The label 
(L) indicates that Acp refers to vector mesons with longitudinal polarization 
only. 



VV 


B{B ^ VV)/10-^ 


h 


B{B ^ VlVl)/W-' 


Acp 


p+p- 




n q7s+o.025 


23.7 ±3.2 


0.06 ±0.13 (L) 




0.68 ±0.27 


0.70 ±0.15 


0.48 ±0.21 


-0.4 ±0.9 (L) 


p-p' 


18.2 ±3.0 


n 012+0-044 


16.6 ±2.9 


-0.08 ±0.13 


P^UJ 


10.6+1^ 


0.82 ±0.11 


8.7 ±2.4 


0.04 ±0.18 




1.28l°:i [25] 


0.80^°:}^ [25] 


1.02 ±0.32 






5.6 ± 1.6 


0.57 ±0.12 


3.2 ± 1.1 


0.09 ±0.19 


K*-p+ 


< 12 




< 12 




K*-p^ 


< 6.1 


96+°-°^ 


< 6.1 




K*°p~ 


9.2 ± 1.5 


0.48 ±0.08 


4.4 ± 1.0 


-0.01 ±0.16 


K*y 


9.5 ±0.8 


0.484 ± 0.034 


4.6 ±0.5 


-0.01 ±0.06 


K*-(f) 


10.0 ± 1.1 


0.50 ±0.05 


5.0 ±0.7 


-0.01 ±0.08 




l.Strr [26] 


O.bGt'oit [26] 


1.0 ±0.7 





Because of similarities in the flavour structure of Tp and Tp"'"' , in all amplitudes 
the coefficient 63 appears together with the factorization coefficient 04 in the combination 
+ Bvj^V2/^ViV2h- This has been noted before in the context of PP and PV final states 
[23]. 

4 Experimental results and input parameters 

Available data on the decays of B~ and B^ mesons into a pair of light vector mesons 
are displayed in Tables 1 and 2. The results are from [24] unless indicated otherwise. 
CP averaging is understood for branching ratios and Jl- No data are available yet on 
Bg — > VV decays. 

Table 3 collects the input parameters used in our analysis. The values of my, Fy, 
'nT'B, \Vus\ and \ Vcb\ have been taken from [28]. They have only small uncertainties, which 
we neglect. Our choice for |146/K&| is compatible with the exclusive determinations 
quoted in [28]. We prefer those over the inclusive values since we use Vub in exclusive 
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Table 2: Experimental results for further B — > VV decays [24]. Quoted are 
the CP-avcraged branching fractions in units of 10~^. The B meson is either 
a B^ or a B^. 







p-0 [27] 




K*+K*- 


< 3.4 


< 71 


< 3 


< 2.7 


< 141 




UJU) 


pV [27] 




<t>4> [27] 


< 1.5 


< 4.0 


< 0.33 


< 1.2 


< 0.2 



processes where the form factors rely on similar theoretical methods (light-cone QCD 
sum rules, lattice) as in the exclusive extraction of \Vub\- The lifetimes of B~ and B^ 
are also from [28]. On the other hand, the lifetime of Bg is put equal to rg^, which 
is expected theoretically to hold to very high accuracy. The value of r^^ from [28] is 
compatible with this, but is still affected by a larger error. 

The number for sin 2/3 is the average of [24] from CP violation in 6 ^ ccs modes. 
The angle 7 corresponds to the result of global CKM fits [29,30]. The standard model 
parameters m6, mc, mt and Mw are the same as in [6]. Changes in these values 

have been small in comparison with the relevant uncertainties. The quark masses are 
running MS-masses. 

The decay constants fv can be determined from data onV ^ and r — > Vu. We 
use the values quoted in [23] . The transverse decay constants fy need to be computed 
theoretically, for instance with QCD sum rules. The results we use for fy have been 
compiled in [31]. The B ^ V form factors are from QCD sum rules on the light cone 
[32]. These results do not yet incorporate some improvements in the treatment of SU{3) 
breaking that has been achieved in the meantime (see comments in sec. 2.3 of [1]). 
The uncertainties on the form factors in Table 3 are taken to be somewhat larger than 
reported in [32]. The Gegenbauer coefficients 0^2 are still rather uncertain. We adopt 
numbers of the typical size found in QCD sum rule calculations [32], [33] and allow 
for sizable uncertainties. The range of numbers for the 5-meson decay constants is 
representative of results from recent unquenched lattice simulations (see sec. 2.4 in [1] 
for a review and detailed references). The parameter Ab is not well known at present. 
We shall consider here the generous range already used in [6]. No attempt is made to 
account for SU (3) breaking in this quantity. 
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Table 3: Input parameters for B VlVl decays. Here Bq stands for either 
B~ or Bd- The values of the scale dependent quantities fy = fyilA given 
for = 1 GeV. The scale dependence of 0(2 is neglected. 



Light vector mesons 


V 


my/MeV 


Ty/MeV 


/y/MeV 


/^/MeV 


aX 


aX 


P 

K* 




776 
783 
894 
1019 


149 
8 
51 

4 


209 
187 
218 
221 


165 ±9 
151 ±9 

185 ± 10 

186 ±9 





0.1 ±0.1 



0.1 ±0.3 
0.1 ±0.3 
0.1 ±0.3 
0.1 ±0.3 


B mesons 


B 


ms/GeV 


Te/ps 


/fi/MeV 


Aij/MeV 


B- 


5.28 


1.64 


200 ± 30 


350 ± 150 


Bd 


5.28 


1.53 


200 ± 30 


350 ± 150 


Bs 


5.37 


1.53 


230 ± 30 


350 ± 150 


Form factors 




Ao"'^-(0) 




<^-^*(0) 




0.30 ±0.04 


0.28 ±0.05 


0.37 ±0.05 


0.36 ±0.05 


0.47 ±0.06 


SM parameters 


Agi/MeV 


mb(mfe)/GeV 


mc(mfe)/GeV 


mt{mt)/GeV 


Mv^/GeV 


225 


4.2 


1.3 ±0.2 


167 


80.4 


\Vus\ 


\Vcb\ 


\Vub/Vcb\ 


7 


sin 2/5 


0.226 


0.0416 


0.09 ±0.01 


(67 ±12)° 


0.681 ±0.025 



5 Phenomenological analysis 
5.1 B — > VlVl branching fractions 

The branching fraction of a decay B — > VilV2l is obtained from the corresponding 
amplitude A as 

B{B ^ VilV2l) = S -p— \A{B ^ V^lV2l)? (120) 

Here is a symmetry factor with 5" = 1/2 if and V2 are identical and S — 1 otherwise. 

Predictions of CP averaged branching ratios are compiled in Table 4 and Table 5 
for strangeness-conserving and strangeness-changing B — > VilV2l decays, respectively. 
Absolute branching fractions have in general sizable uncertainties from hadronic input 
quantities, for instance from B form factors. Taking ratios or other combinations 
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Table 4: CP-avcragcd branching fractions for B VlVl decays with AS* = 0. 
The sensitivity to variations in the input parameters according to Table 3 is 
displayed where the upper (lower) entry corresponds to the larger (smaller) 
value of the parameter. The renormalization scale /i is varied between 2mb 
and m(,/2. The model parameters Xa,h{pa,h-i4>a,h) from power corrections 
are varied within the range given by < pa,h < 1 and < 0a,_h" < Stt. Here 
upper (lower) entries refer to positive (negative) lm.XA,H- The appropriate 
units for each mode are given in square brackets. 



mode 


central 




aX 


Ab 


Ib 


P 


Xa 


Xh 


V u 
y CO 


7 


B p-p^ 


17.5(10-6 




+4.5 
-4.0 


+1.9 
-1.4 


-0.9 
+2.5 


+0.5 
-0.5 


+0.1 
-0.0 





+1.5 
-1.4 


+4.1 
-3.6 


-0.2 
+0.2 


B^ — i> p~u 


15.5(10-6 




+5.6 
-4.7 


+1.4 
-0.9 


-0.8 
+2.0 


+0.3 
-0.3 


+0.3 
-0.3 


-1.2 
+0.9 


+1.2 
-1.1 


+3.4 
-3.0 


-0.8 
+0.7 


B- p-ct> 


6.0(10-9] 


+0.9 
-0.8 


+3.7 
-2.2 


-1.5 
+4.8 


+0.9 
-0.8 


-1.1 
+3.4 




+2.8 
-2.2 


+0.0 
-0.0 


+1.1 
-1.0 


B- K*-K*^ 


2.7(10-^ 


+0.9 
-0.8 


-0.8 
+0.9 


+0.2 
-0.4 


-0.1 
+0.1 


-0.4 
+0.5 


-2.5 
+3.9 


-0.3 
+0.3 


+0.0 
-0.0 


+0.5 
-0.5 




3.3(10-^; 




+0.3 
-0.3 


+4.3 
-1.6 


-1.3 
+6.1 


+0.7 
-0.6 


-0.0 
+0.7 


+2.2 
-1.7 


-2.2 
+3.7 


+0.8 
-0.7 


+0.3 
-0.3 


Bd ^ p^cc; 


8.0(10-8; 




+3.1 
-2.5 


-3.3 
+3.9 


+1.1 
-2.4 


-0.4 
+0.5 


-1.4 
+1.4 


-2.7 
+22.5 


-1.6 
+1.8 


+0.4 
-0.3 


+2.3 
-2.1 


Bfi — > ujuj 


5.0(10-^; 




+0.7 
-0.6 


+3.7 
-1.7 


-1.4 
+5.4 


+1.0 
-0.9 


-0.2 
+1.0 


-1.3 
+2.1 


-2.3 
+3.3 


+0.9 
-0.8 


-0.4 
+0.4 


Bd pV 


2.8(10-9; 




+0.4 
-0.4 


+1.7 
-1.0 


-0.7 
+2.2 


+0.4 
-0.4 


-0.5 
+1.6 




-1.0 
+1.3 


+0.0 
-0.0 


+0.5 
-0.5 


Bd uj(() 


2.4(10-9] 


+0.5 
-0.4 


+1.4 
-0.9 


-0.6 
+1.8 


+0.3 
-0.3 


-0.4 
+1.2 




-0.9 
+1.1 


+0.0 
-0.0 


+0.4 
-0.4 


Bd ^ p+p- 


25.8(10-6 




+7.6 
-6.6 


-2.0 
+1.6 


+1.0 
-2.4 


-0.3 
+0.3 


-0.3 
+0.1 


+2.4 
-1.8 


-1.5 
+1.5 


+5.8 
-5.2 


-0.9 
+0.8 


Bd K*^K*^ 


3.2(10-^ 


+0.9 
-0.8 


-0.9 
+0.9 


+0.2 
-0.4 


+0.0 
-0.0 


-0.7 
+0.9 


-2.5 
+3.0 


-0.3 
+0.4 


+0.0 
-0.0 


+0.5 
-0.4 


Bs K*^p^ 


5.6(10-^; 




+0.4 
-0.4 


+6.7 
-2.7 


-2.1 
+8.8 


+1.2 
-1.0 


-0.0 
+1.1 


-0.4 
+1.5 


-4.2 
+7.5 


+1.3 
-1.2 


+0.4 
-0.4 


Bs K*^oj 


6.5(10-^; 




+0.9 
-0.7 


+5.5 
-2.3 


-1.9 
+7.6 


+1.0 
-0.9 


-0.0 
+0.9 


-0.9 
+1.6 


-3.9 
+6.6 


+1.2 
-1.1 


-0.5 
+0.4 


Bs K*^(t) 


3.4(10-^] 


+1.1 
-1.0 


-1.3 
+1.4 


+0.3 
-0.7 


-0.2 
+0.2 


-0.6 
+0.6 


-3.3 
+5.1 


-0.7 
+0.8 


-0.0 
+0.0 


+0.5 
-0.5 


Bs ^ K*+p- 


37.2(10-6 




+11.8 
-10.2 


-2.9 
+2.3 


+1.4 
-3.3 


-0.6 
+0.6 


+0.0 
-0.3 


-0.2 
+0.4 


-2.6 
+2.7 


+8.4 
-7.5 


-1.2 
+1.1 


Bd K*+K*- 


2.9(10-8] 








+0.9 
-0.8 


-1.0 
+2.2 


-2.8 
+19.1 




+0.5 
-0.4 


-0.4 
+0.3 


Bd^(t>(t> 


2.5[io-9; 










+0.8 
-0.7 


-1.3 
+3.4 


+16.7 
-2.4 




+0.0 
-0.0 


+0.5 
-0.4 
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Table 5: CP-averaged branching fractions for B — > VlVl decays with AS — 1. 
The sensitivity to variations in the input parameters according to Table 3 is 
displayed where the upper (lower) entry corresponds to the larger (smaller) 
value of the parameter. The renormalization scale /i is varied between 2m(, 
and mb/2. The model parameters Xa,h{pa,h, (Pa,h) from power corrections 
are varied within the range given by < pa,h < 1 and < (f)A,H < 27r. Here 
upper (lower) entries refer to positive (negative) lmXA,H- The appropriate 
units for each mode are given in square brackets. 



rnnnp 


central 






"2 


Ad 
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of suitable branching fractions can eliminate part of the uncertainties and lead to theo- 
retically cleaner observables. In spite of this it is still interesting to present the theory 
expectations for the branching fractions, which can be directly confronted with exper- 
imental data. In addition, wc use Tables 4 and 5 to display in detail the sensitivity of 
the results on the most important input parameters. 

In Fig. 5 we compare theory and experiment for B — > VlVl branching fractions, for 
which measurements are available. For this comparison the form factors have been fixed 
to their central values. In the present discussion we will assume that the Standard Model 
is valid. Under this assumption the comparison with experimental data will serve as a 
test of the theory of QCD interactions in hadronic weak decays. It should be kept in 
mind that possible deviations between predictions and measurements may in principle 
indicate the existence of New Physics. In order to disentangle New Physics from QCD 
effects it is important to consider observables with very little hadronic uncertainty. We 
will discuss several examples for this in the following sections. For the moment we stay 
with the discussion of the theoretically less clean absolute branching fractions for the 
purpose of testing the method of QCD factorization, under the (provisional) assumption 
that physics beyond the Standard Model is absent. 

Fig. 5 shows good agreement of theory and experiment within errors. An exception 
is B~ — > p~iJ, where the measured branching ratio is somewhat low with respect to the 
expectation from theory. The reason could be an overestimate of the B ^ uj form factor, 
or a statistical fluctuation. The only other measured decay with an a; is — > K*^uj. 
Here the theory result is also above the central experimental value, but the uncertainties 
in the latter are still larger than for B~ — > p~uj and prevent any firm conclusion. 

We emphasize that theory and experiment agree very well in the three p-meson 
channels Bd — >■ p'^p~, B^ p^p^ and B~ P~ P^ ^ as it has also been found in [11]. 
The channel is a colour-suppressed mode and comes with large uncertainties. Hard 
spectator scattering plays an important role and therefore the sensitivity to the poorly 
known parameter is large. Still the experimental result can be accounted for naturally 
with default values of the hadronic parameters. It is remarkable that the observed pattern 
of all three pp branching ratios, which exhibit rather different values, is nicely reproduced 
within QCD factorization. The fact that this works for the central B ^ p form factor 
Aq supports the numerical value used for this quantity. 

The penguin modes B4 x*^K*^ and B^, K*^p^ tend to have relatively smaU 
predicted branching ratios, which however stretch into the range of measured values 
within errors. The compatibility is better for B^ ^ ^*oj^*o ^j-^g^j^ fQj^ _^ K*^p^ . 
At the same time the latter mode is also seen to be very sensitive to the annihilation 
contributions. 

It is interesting to note that the central values of the experimental and theoretical 
results are particularly close for the penguin decays B~ K*^p~ , Bd — > B~ 
K*~(j). On the other hand, the dependence on weak annihilation is very strong. The huge 
variations from these effects shown in Fig. 5 suggest that, at least for these channels, 
the annihilation model used by us is likely to overestimate the related uncertainty. 

Further branching ratio predictions and information on the various error sources for 
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Figure 5: Comparison between theory predictions [dots (black), left bar (cyan), middle 
bar (marine blue)] and experimental results [right bar (orange)] for Bd VlVl modes, 
for which measurements are available. The theoretical error bars display the hadronic 
errors without [left (cyan)] and with [middle (marine blue)] the model-dependent error 
estimate for annihilation topologies. The form-factor uncertainties are not included in 
the error bars. The black dots are the central values of the theory predictions where all 
power corrections have been omitted. From experiment only upper limits are known for 
the two K*~p channels. 
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all 34 5 — > VlVl decays can be obtained from Tables 4 and 5. 

Our results include estimates of some effects that are suppressed by a factor of 
^QCD/iTib- These corrections are weak annihilation and the effects proportional to 
(see eq. (25)). Terms at this order arc not calculable in QCD factorization. They have 
still been included as model estimates in order to permit us to assess the sensitivity of 
factorization predictions on potentially important power corrections. Weak annihilation 
is the most prominent example. For the default choice of input parameters the impact of 
power corrections on the predicted branching ratios is in general small. This can already 
be seen from Fig. 5, where central results with all power corrections omitted are indi- 
cated by the black dots. They differ very little from the central theory predictions that 
include such effects. To make these statements more quantitative, we list the differences 
between the central values for all AS — branching ratios without and including power 
corrections, BR{n.o power corr.)/Si?(default) — 1, in the order of appearance in Table 
4, in %: 

-1, +2, -6, +6, +15, +7, -27, -6, -6, -5, -16, -12, -9, +6, +1, -100, -100 

(121) 

The same information for the A^" = 1 decays of Table 5 reads 

+5, +5, +3, +7, +9, +12, +6, +7, +2, +0, -8, -17, -19, -100, -100, -100, -100 

(122) 

The deviation is —100% for the six pure annihilation decays, which have no leading- 
power contribution. In all other cases the impact of the default power corrections is 
rather moderate or indeed very small, notably for the dominant decay channels. 

We finally comment on the impact of the long-distance electromagnetic penguin cor- 
rection defined in (39) and discussed in appendix A. This contribution affects only decays 
with the emission of ou or (j), where it enters through the coefficient a" + Og. The 
long-distance effects are sizable, on the scale of this coefficient, for p° and u, but much 
less in the case of 0. Since the long-distance terms are of order a = 1/129 their overall 
contribution is in general very small. This is particularly true for the AS — 1 decays 
where the up-quark sector is also CKM suppressed. For the AS — transitions the 
absence of the term in (39) would change branching ratios at the level of a few percent 
at most and below the size of most of the other uncertainties. The situation is similar 
for the direct CP asymmetries in the A^* = modes with the exception of p~p^, p~4', 
p'^u, p^(j) and tc;0, where the impact is relatively large. However, in any case, the direct 
CP asymmetry is very uncertain for p'^uj and it is very small for the remaining channels. 

5.2 Direct CP violation in B ^ VlVl 

Direct CP asymmetries require the presence of a strong as well as a weak phase difference 
between two interfering amplitudes. In the heavy-quark limit this phase difference arises 
at order a^. It is therefore parametrically suppressed and at the same time sensitive 
to uncalculable power corrections. This makes it difficult to obtain accurate predictions 
for direct CP violation. At present the most precisely measured direct CP asymmetry 
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in B decays is Acp{B — >■ K~^7r~) = — 0.097 ± 0.012 [24]. The corresponding strong 
phase difference is small (~ 15°) [34], but has a sign opposite to the 0{as) result in the 
heavy-quark limit. This may indicate the importance of A/ mi, corrections. 

In Tables 6 and 7 we show estimates of direct CP asymmetries for the decays under 
discussion. The values have large uncertainties, as anticipated. Most of the asymmetries 

Table 6: CP asymmetries for B — VlVl decays with AS = 0, defined as 
Acp = {B{B ^ /) - B{B ^ f))/{B{B -> /) + B{B ^ /)). The sensitivity 
to variations in the input parameters according to Table 3 is displayed where 
the upper (lower) entry corresponds to the larger (smaller) value of the pa- 
rameter. The renormalization scale is varied between 2mh and mh/2. The 
model parameters Xa,h{pa,h,4>a,h) from power corrections are varied within 
the range given by < pa,h < 1 and < (I)a,h < 27r. The appropriate units 
for each mode are given in square brackets. We refrain from quoting estimates 
of CP asymmetries for pure annihilation modes. 
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are small or moderate, but there can be exceptions. Large asymmetries may occur when 
the interfering amplitudes have comparable magnitude and a substantial strong relative 
phase. Examples are the AS* = decays with a colour suppressed tree contribution 

(~ 02), as Bd — > p^p^, ujuj or Bg — ^ K*^p^, K*^uj. Despite the as factor the strong phase 
difference can here be naturally more sizeable. Generically, a decay amplitude of the 
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Table 7: CP asymmetries for B — > VlVl decays with AS — 1 (see caption of 
Table 6 for details). 
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form 

A{B Ml Ma) ~ e-'^ - pe^"^ (123) 
with p, (p, 7 real, leads to the direct CP asymmetry 

^ 2psin0 sin 7 

1 + — 2p cos 4> cos 7 

For sin 7 ^ 0.92, a value p = 0{1) and a substantial phase give a large asymmetry. In 
the case of B^, PlPl the central values p — 0.36, cf) — 49° give Acp — 53%. 

5.3 Sensitivity to a;-0 mixing 

In the other sections of this paper the vector mesons (f) and u arc always implemented 
as pure ss and {uu + (iJ)/^/2-states, respectively. Here we investigate the sensitivity of 
our results to the deviation from this case of ideal mixing. We assume that other effects 
with Zweig-rule suppression are negligibly small. We neglect, for example, Zweig-rule 
forbidden matrix elements of the type {(l){ss)\{ub)v-A\B'~) . 



(124) 
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Mixing can be introduced by the following parametrization: 



0(1020) ^ss cos g + ""^^^^ sing (125) 

V2 

a;(782) ^ '^''^ "^"^ cos 6 - ss sin 6. (126) 
v2 

The ideal mixing angle in this parametrization is ^ = 0. According to sum-rules quadratic 
in meson masses [28], the mixing angle can be estimated to be ^ = 3.4°. The results 
of varying the mixing angle up to ^ = 6.8° are shown in Table 8. For most branching 



Table 8: Dependence of i? — VlVl branching fractions on a;-0 mixing. The 
variation of the branching fractions is given for two values of the mixing angle 
9. The upper (lower) value corresponds to = 6.8° {9 = 3.4°). 



mode 


default value 


deviation 


mode 


default value 


deviation 


B~ — > p~U! 


15.5[10-s 




-0.2 
-0.0 


B- K*-uj 


1.7[10-6 




-0.4 
-0.2 




6.0[10-9] 


+207.5 

+49.8 


B- K*-(t) 


4.1[10-6 




+0.4 
+0.2 




8.0[10-''] 


+0.2 
+0,1 


D,, A" '"^' 


i.i[io-«] 


-0.,'') 
-0.:i 




5.0[10-'] 


-O.i 
-0,0 




3.7[10-'] 


+0.,^) 
+0.:-! 




2. 8 [10-"] 


-1,!) 
-1.2 


D„ fPo 


5. 9 [10-'] 


-0.1 
-0.0 


Bd uj(f) 


2.4[10-9 




+10.0 
+1.7 


Bg — > uj4> 


4.4[10-^] 


+28.7 
+4.6 


B, K*^uj 


6.5[10-^ 




-0.6 
-0.3 


Bs^(t>(t> 


15.5[10-6 




-0.3 
-0.0 


Bs K*^(t) 


3.4[10-^ 




+0.6 
+0.3 


Bg — > UJUJ 


3.2[10-«] 


+1.6 
+0.3 




2.5[10-'' 




-0.4 
-0.3 


Bg ^ pOu; 


1.5[10-9 




+11.1 
+3.5 



fractions the effect of a nonvanishing mixing angle 9 ~ 3.4° is very small, in particular for 
the important modes B~ — > p~uj, B~ K*^(j), Bg (jxp. On the other hand, the modes 
Bfi p^(f), Bfi — > Bg — > Bg — > p^cu have a significant dependence on deviations 
from ideal mixing. The largest effect is observed for B~ — > In this case p~uj 
feeds into the former channel through mixing with a more than three orders of magnitude 
higher branching ratio compared to B~ p~(j){ss), which overcompensates the small 
mixing angle: 

B{B- p~0)mix ~ sin^ 9 B{B- p'u;) (127) 

A recent discussion of hadronic B decays, mostly with charm in the final state, for which 
mixing has a large impact can be found in [35]. Their estimate of B" — > is 
compatible with ours. 
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5.4 Unitarity triangle from CP violation in Bd — ^ PlPl 
5.4.1 Determination of p, 77, 7 and ct 

The time dependent CP asymmetry in Ba — > is given by 

(128) 

The parameters Sp and Cp have been measured to be 

5p = -0.05 ± 0.17 Cp = -0.06 ±0.13 (129) 

as quoted by [24], based on results of BaBar [36] and Belle [37]. Together with the 
experimentally well determined quantity sin 2/3 from CP violation in B ipK^ decays, 
the parameter Sp can be used to fix the CKM unitarity triangle. The value of sin 2/3 
from Table 3 implies /3 = (21.5 ± 1.0)° or 

r = cot/3 = 2.54 ±0.13 (130) 

In terms of the improved Wolfenstein parameters p and 77 [38] the unitarity triangle is 
then determined by 

p^l-Tfj (131) 



1 

77 = 



1 + t5'p)(1 + rpCOS0p) 



1 + T^)S, 

-^(1 - 5-2) (1 + cos 0^)2 - 5-^(1 + T^){Sp + sin2/3)r2sin 



2 cm ^ 



(132) 



These formulas have been derived in [19,20] for B tt+tt^, but they apply to the case of 
B P~^P2 as well. The parameters Vp and 4>p are hadronic quantities. They are defined 
here through 

a, + at + a«o + (h + h + 2h - + i^D 
where all coefficients a^, 6j refer to the p^p^ final state and 

-^-|^ = ^#^-5.10- (134) 
^pp m^^o (0) 

The real quantities and 0p are the magnitude and phase of the penguin-to-tree am- 
plitude ratio in B ^ PlPl- They are independent of CKM parameters. Numerically we 
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find 



Vp = 0.038 ± 0.005 (/X, a^) 1°;°^^ (p^, (135) 

0p = 0.23 ± 0.09 (me, a^) 1^;?^ (pa, 0a) (136) 

cos cj>p = 0.037 ± 0.005 (//, «^) 1^26 (P^, 0a) (137) 

The first error is from the uncertainties in the input parameters ^40, q;2) fpi ^b, Jb, 
rric and a variation of the renormahzation scale /i between m6/2 and 2m6 around its 
default value /i — mb. The dominant sources of uncertainty are indicated in brackets. 
The second error reflects the sensitivity to the parameters p^, 0a, Ph and (()h used 
to model power corrections from weak annihilation (A) and in the spectator scattering 
amplitude (H). We have used < pa,h < 1, < 4>a,h < 27r. The second error is entirely 
determined by weak annihilation. 

The phase 0p is parametrically suppressed since it arises only at order ag or Aqcn/'iTT'b- 
Its precise value is rather uncertain, in particular due to the model dependence of power 
corrections, which may compete numerically with the calculable 0{as) term. Fortunately 
the dependence of fj in (132) on 0p is very weak [19,20]. In addition, Vp is a small 
parameter, even smaller than the corresponding quantity r„ in Ba tt'^tt' . The smaller 
size of the penguin contribution in the case of vector mesons as compared to pseudoscalars 
has been pointed oTit before in the context of QCD factorization [9,23]. The formulation 
in (132) makes it particularly transparent to analyze the impact of a small penguin 
correction on the determination of the unitarity triangle. To linear order in r^, eq. (132) 
implies the simple relation 

1 + tS — /I — 5*2 
^ = (l+r^)Sp ' + ""^ ''''' '^^^ ^^^^^ 

In this approximation fj and p depend only on the real part of the penguin-to-trec ratio. 
As can be seen from (132), second order corrections in Vp are further suppressed by 
sin^ (pp. 

With fj and p also the CKM angles 7 and a can be computed: 

fj 

7 = arctan -, a = tt — /? — 7 (139) 

1 — TT] 

It is instructive to write down the expressions for small values of Sp, which are suggested 
by the data in (129). To first order in both Sp and we find 

7 = arctan t + + rrp cos 0p (140) 

a = I - Y -^»^pcos0p (141) 

For S'p = and in the absence of a penguin contribution one has a. = 90° and 7 = 
(68.5 ± 1.0)°. Non-zero values of the observable Sp and the theoretical quantity rpCOS0p 
then compete in shifting 7 and a away from these lowest-order approximations. 
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Evaluation of the exact formulas (132) and (139) gives 

fj = 0.350 ± 0.013 (r) ± 0.012 (Sp) ± 0.008 (r^ cos 0^) (142) 

7 = 72.4°± 1.3°(r) ±5.1° (Sp) ± 3.2° (r^ cos 0^) (143) 

Nearly identical results are obtained for 7 when the first order expression (140) is em- 
ployed. The approximations (140) and (141) work to very good accuracy in the relevant 
range of Sp and r^. This greatly facilitates the determination of 7 and a and the analysis 
of errors, which can simply be read off from (140) and (141). 

The calculation of 7 in [11] using the longitudinal part of the time dependent CP- 
asymmetry in the p"'"p~-system and /3 as input yields a similar result for the hadronic 
error of ±3°. 

The determination of 7 in (143) is considerably more precise at present than mea- 
surements using B — > DK tree-level decays. Belle has found [39] 

7 = {53tll{stat) ± 3{sys) ± %model))° (144) 

and a recent analysis from BaBar gives [40] 

7 = (76 ± 22{stat) ± 5{sys) ± b{model))° (145) 

Within the errors both are well compatible with (143). 



5.4.2 Bounds on UT parameters 

Useful information on the angle 7 can also be obtained in the form of a lower bound, 
which is even less sensitive to theory input than the result in (143). It relies only on 
the conservative condition that rcos0 > 0, which holds in the heavy-quark limit. This 
bound has been derived in [19,20]. Further discussions may also be found in [41,42]. The 
bound is vahd as long as > — sin 2/3 and reads 



7 > arctan ^ ' (146) 

The constraint (146) can be evaluated using CP violation in B vr+vr^ {S = Sj^) or 
in B ^ pXpJ^ {S = Sp). The derivation of (146) is identical for both cases. In fact, 
since the (positive) penguin correction r^cos^p is smaller than r^rCos^^, the bound is 
expected to be more stringent using Sp instead of S^,-. This expectation is indeed bourne 
out by the experimental result Sp > S-„ = —0.61 ±0.08 (see sec. 5.7), which implies that 
Sp gives the better constraint. Qualitatively, these features can also be understood from 
the approximate relation (140). 

To linear order in S the bound (146) becomes 

g 

7 > arctan r+— (147) 
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in agreement with (140). 

Using S = Sp ^ -0.05 (central), -0.22 (la), -0.39 (2a), we obtain from (146), 
respectively 

7>67°, 62°, 57° (148) 

The linear approximation (147) gives practically identical results. We remark that the 
relevant values of Sp fulfill the condition Sp > — sin 2/3, under which the bound can be 
applied. 

The penguin correction is expected to shift the numbers in (148) by approximately 
+6° to yield the actual value of 7. The bound is therefore quite stringent. Within 
a Standard Model interpretation it eliminates already a sizable fraction of the allowed 
range from the direct measurements in (144) and (145). 

Bounds similar to the one for 7 can also be derived for fj and p [19,20]. The lower 
bound for fj is given by the right-hand side of (138) with Vp put to zero, the upper bound 
on p then follows from p — 1 — rfj. With the same input for Sp as in (148) we find 

f^> 0.338, 0.326, 0.314 (149) 

and 

p< 0.143, 0.172, 0.203 (150) 



5.4.3 Precision determination of \Vub\ from sin 2/3 and Sp 

The preceding analysis has a further interesting application regarding the determination 
of \Vub\ from sin 2/3 and Sp. The value of \Vub\ determined in this way may be affected 
by New Physics entering CP violation in ipKs and PlPl- "^^^ presence 

of non-standard contributions can be revealed by comparing the extracted value of \Vub\ 
with the result for jV^^I from an independent method. An important example is the 
direct determination of \Vub\ from semileptonic, exclusive or inclusive, b — > ulu decays, 
which are most likely independent of physics beyond the Standard Model. It is clear 
that the usefulness of such a New Physics test will depend on how precisely \Vub\ can 
be determined. We will show that sin 2/3 and Sp offer a particularly clean and accurate 
determination of \Vub\- 

The magnitude of \Vub\ is proportional to Rb = \/ p'^ + ff- Using the exact formulas 
in (131) and (132), we expand in Sp and Vp. This is motivated by the smallness of 
the theoretical parameter Vp and the empirical observation that also Sp is small, as we 
have discussed in section 5.4.1. Treating Sp and r-p as small quantities of the same order 
we find 



Rb^VW+W-^r=^ (151) 

Vl + T"' 

1 + ^ Y + rrp COS + ^ ^ Y + TTp cos {Sp cos 0^ + 2Trp sin^ 0p) 
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where we have neglected terms of the fourth order. Through terms of third order in Sp 
and rp (5^, S^^r-p, Sprj, r^) eq. (151) is exact. 

The basic features of (151) are easy to understand from the geometry of the unitarity 
triangle. If 5"^ = Tp = then a — tt/2. In this case Rb — sin/3 = l/Vl + r^, which gives 
the leading term in (151). Because sin [3 is the minimum value that i?5 can take for fixed 
(3, first order corrections in Vp and Sp are absent and the second-order term is strictly 
positive. The protection of (151) from first-order corrections in Sp and is the basis for 
a precise determination of Kb- 

The quantity Sp/2 + rvp cos (pp appeared already in (140), (141). For Sp < there 
is a further cancellation in this term with the penguin shift rVp cos (f)p. Taking Sp = 
-0.05±0.17 (129), T = 2.54±0.13 (130) and the conservative range cos0p = 0.04±0.03 
we have 

^ + rrp cos 0p = 0.077 ±0.114 (152) 

The range of Vp cos (pp covers the result obtained from the QCD analysis in sec. 5.4.1. As 
we will show in sec. 5.6, r^cos^p can also be determined by independent experimental 
information on the penguin mode — > K^^K^, which confirms the values employed 
here. 

Through second order in Sp and the correction factor relative to the lowest-order 
result in (151) reads, using (152), 



i + i(^ 

2 V 2 



1 003+°°^^ 

l.UUO_o,oo3 



(153) 



We remark that the lower hmit of 1 for this factor is an absolute bound. The third-order 
term in (151) is less than about 0.2rpSp/2 ^ 0.2 • 0.04 • 0.1 ^ 0.001 and thus completely 
negligible. 

Using sin/3 = 0.366 ± 0.016 from Table 3 we obtain 



Rb — sin (3 



. 1 fSp 

1+2 ( Y + ^'^pCos, 



0.367 ± 0.0161^02 



(154) 



From [28] we have 

A= |V;.| = 0.226 ±0.002 

This imphes 



iKbl = 0.0416 ± 0.0006 



(155) 



IK 



ub\ 



VoMcd 



ud 



Rh 



A 



1 - 



A2 



Rb iKbl = (SMt^o fM ± 0.05(K6) ± 0.03(K.)) • 10 



-3 



(156) 

The uncertainty is dominated by the error in P ((±0.15) • 10~^), followed by the error in 
the correction from Sp, rp ((l^oD ' 10~^) and the error in Kb- Adding errors in quadrature 
the final result reads 

iKbl = (3.54 ±0.17) - 10-3 (157) 
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It corresponds to a ratio iKib/Kbl = 0.085 ± 0.004, in agreement with Table 3. The 
value in (157) should be compared with the direct measurements of \Vuh\ in 6 — > ulu 
transitions. A recent analysis of inclusive decays gives [43] 

\Vub\ = (3.70 ± 0.32) • 10-3 (158) 

The exclusive determination from B — > nlu decays has been investigated in [44] with the 
result 

iKftl = (3.36 ±0.23) ■ 10-3 (159) 

Related discussions, in the context of QCD sum rules, can be found for instance in [45] 
and [46]. Using an average of data from lattice QCD, [47] quotes for the determination 
from exclusive decays 

= (3.54 ± 0.40) • 10-3 (160) 

The results in (157), (158), (159) and (160) are in very good agreement with each 
other. They provide us with a test of Standard Model CP violation in S — > t/jKs and 
B — > p'^p~ (157) against the \Vub\ determination from tree-level, semileptonic b — > ulu 
decays (158), (159) and (160). 

Numbers for \Vub\ very similar to (157) have been obtained from global fits of the 
unitarity triangle performed by the CKMfitter [29] and UTfit [30] collaborations, which 
quote 

\Vub\ = (3.57 ± 0.17) • 10-3 (CKMfitter) (161) 
\Vub\ = (3.55 ± 0.15) • 10-3 (UTfit) (162) 

While such global fit results summarize our current overall knowledge of quark-mixing 
parameters, they do not exhibit explicitly the individual pieces of information that deter- 
mine this knowledge. We emphasize here that the precise result in (157) can be obtained 
from sin 2/5 and Sp alone, with only very moderate requirements on the accuracy of the 
penguin contribution ~ Vp from theory. The representation proposed in (154) makes this 
statement particularly transparent. 

The result in (157) is currently the most precise determination of \Vub\- Since the error 
is dominated by the uncertainty in sin f3, an even higher precision will be achieved by a 
more accurate measurement of sin (3 as it is expected at the upcoming LHC experiments. 
For instance, with a determination of sin/3 to 1%, the error in (157) would shrink to 
±0.08 • 10-3, corresponding to a precision of 2% for \Vub\- 

Constraint on New Physics phase in Ba — mixing 

The preceding analyses rely on the assumption of a Standard Model phase in B^ — B4 
mixing. We would hke to examine the effect of a small New Physics phase entering 
only in Bd — -B^-meson mixing [48,49]. In this scenario the New Physics phase shall 
not violate unitarity of the Standard Model CKM matrix. The modified mixing phase 
j3 + Vl, with the New Physics contribution Vl, enters the analysis in the determination 
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of r 



= cot(/5 + fl) from Bd J/ipKs and through mixing-induced CP violation in 
p^p^ . The relation (131) for p will depend now on Q: 



(1 — TTj) — {t + 1]) tan Q 
1 — r tan Q 



(163) 



The measurement of Sp determines r) up to the mixing phase fl. The new relation for rj 
in an expansion in Vp and Sp reads 



V 



r(l-rtanfi) ^ (1-TtanQ) (^rr^ cos 0^ + ^(1 - r tan Q)) 



1 + t2 



l + r2 



(164) 

We note that Q enters the leading term in the expansion with an enhancement of r ^ 2.5. 
The first order term is suppressed relative to the leading order term by a factor of ~ 0.1. 
Using the exact relations for p and fj one finds the following expansion for Rj,: 



Rh 



|1 — rtanf2| 



l-tann(.^I^^ + 



rtanO 



+ 0{rlrpSp,S') 



(165) 



Again we have an enhancement of the dependence on Q by r. For Ri, the error from Sp 
'p is much less important than the error from r (130), in contrast to the case 



and Tp cos 



of (142). 

Relating R}, to \Vub\ as in (156), taking \Vuh\, Sp, r from experiment and Tp and (f)p 
from QCD factorization, the angle Q can be extracted. In general the solution is not 
unique. Here we assume that the new phase Q is small, neglecting discrete ambiguities. 
Such ambiguities may be eliminated with additional measurements. In particular, a 
second solution with large ^2 would imply a negative sign of cos(2(/5 + f^)), which is 
disfavoured by experiment [50]. A more general discussion on the New Physics aspects 
of this analysis can be found in [49] . 

If one disregards solutions with \Q\ > arctan(l/r) ~ 20° and uses the exclusive 
determination (160) of \Vub\, then Q can be determined with an accuracy of few degrees. 
Using the exact expression for Rb{T, Sp, Vp, <pp, fl) we find 



^ - (o.o;?:^(r) triis,) l^:?(r,) lgJ(K.) ;l^(v;.))° 



(166) 



One may note the very small impact of Sp — — 0.05±0.17 and = 0.04±0.03. A different 
value for \Vub\ from other direct determinations such as (158) or (159) leads to very similar 
results. Combining the errors in (166) in quadrature one finds Q = (0.0 ± 2.7)°. 



5.5 Extracting Vp from B K^pj^ 

The precision of CKM angles extracted from CP violation in 5 — > p^pl is ultimately 
limited by our knowledge of the penguin parameters Vp and, to a lesser extent, 0p. Since 
Tp is small, a very moderate accuracy in this quantity is sufficient to obtain a small 
theoretical error for CKM parameters. In [51] it has been proposed to constrain the 
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penguin parameter Vp using the penguin dominated decay B~ — > Kf^p^. We will discuss 
this method in the context of our analysis, comment on the benefits and limitations, 
present an updated numerical evaluation, and compare with the theory results of sec. 
5.4. 

The main idea of [51] is to determine the penguin amplitude from the pure-penguin 
process B~ — > K^p'^ through 

B{B- Kfp,) = ^^11^ \ac{pK*)\' (167) 

Here we defined ap{pK*) = ap{B~ — > K*j^pl) as the coefficient of {iGp/ \^)Xp in the 
amplitude for B" — > K^pJ^, eqs. (67) and (106). They correspond to the charm- and 
up-quark penguin amplitudes for this process. Since au{pK*) and a^pK*) are of compa- 
rable size, and the up-quark amplitude is strongly CKM suppressed, the charm penguin 
completely dominates the branching ratio (167). The penguin amplitude adpK*) can 
be related to the (similarly normalized) penguin amplitude adp) in B^ — * P^P by 
introducing the factor 

acipK*) 



ac(p) 

To lowest order (in and K/mh) this factor would be given by k = f^/fp = 1-04. 
Including QCD corrections this value is reduced to \k,\ = 1.01, and further to \k.\ = 0.97 
by the effects of electroweak penguins. The estimate in (168) includes also the weak 
annihilation terms with default model parameters. Annihilation contributions are thus 
seen to be potentially important. These observations agree with the discussion in [51]. 
In that paper the ratio of the penguin amplitudes in B~ Kf^pJ^ and B^ p'^p~ has 
been parametrized in terms of a factor F, which is related to k through = \/F fK*/fp- 
In [51] a rather wide range for F is assumed, 0.3 < F < 1.5. We will use the same range, 
which corresponds to \k,\ — 0.93 ± 0.36. 

For a given value of [«;[ the penguin parameters r^, (pp are then constrained by the 
ratio 



B{B^ptp^) tb. 



Vr 



cd 



p2 _|_ ^2 _|_ ^2 _|_ 2p fp COS (f)p 



(169) 



where CP averaged rates arc understood. Using r = cot /3, Sp, Cp and the ratio of 
branching fractions in (169) as experimental inputs, the four quantities p, f^, Vp and (pp 
can be determined as functions of A discrete ambiguity in the sign of cos0p can 
be resolved using the heavy-quark limit. The suppression of 0p in this limit singles out 
the solution with cos (f)p > 0. A similar use of the qualitative result cos0p > from 
factorization has been made in [19,20]. Further details on the extraction of p, f], Vp, 
4>p are discussed in 5.6 in the context of a similar analysis. The results for the present 
method are collected in Table 9. Combining and symmetrizing errors we obtain from 
Table 9 

7= (72.5 ±6.9)° a = TT - /3 - 7 = (86.0 ± 7.0)° (170) 
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Table 9: CKM and penguin parameters extracted from r = cot (3 = 2.54±0.13, 
Sp = -0.05 ± 0.17, Cp = -0.06 d= 0.13 and b = B{B- Kfpr^)/B{B 
PlPl) ~ 0.186 ± 0.049. The penguin correction factor is taken to be = 
0.93 ±0.36. 





central 


r 


Sp 




b 




p 


0.110 


-0.010 
+0.012 


-0.030 
+0.031 


-0.002 
+0.028 


-0.005 
+0.005 


+0.010 
-0.024 


V 


0.350 


-0.013 
+0.014 


+0.012 
-0.012 


+0.001 
-0.011 


+0.002 
-0.002 


-0.004 
+0.009 


7[deg] 


72.5 


+1.0 
-1.1 


+5.0 
-5.1 


+0.3 
-4.8 


+0.8 
-0.9 


-1.7 
+3.9 




0.040 


-0.002 
+0.002 


+0.000 
-0.000 


+0.000 
-0.001 


+0.005 
-0.006 


-0.011 
+0.027 




-0.32 


+0.00 
-0.00 


+0.01 
-0.01 


+0.69 
-1.17 


+0.04 
-0.06 


-0.13 
+0.13 



where Sp is the largest source of uncertainty. The results agree very well with those in 
(143). Eq. (170) is an update of the results quoted in [51]. We have checked that we 
obtain the numbers given in that paper if we use the same input. 

A disadvantage of the method just described is that the charm-penguin amplitudes 
in B~ — > and in B^i PlPl ^'^^ related in full QCD by SU{3) flavour 

symmetry alone. The SU{3) argument relating ac{pK*) and ac(p) strictly holds only 
to leading order in the heavy-quark limit. At the level of power corrections from weak 
annihilation these penguin amplitudes are not related by SU{3). This can be seen from 
eqs. (95) and (106), which show that the QCD annihilation penguins are determined 
by the coefficient 63 for ac{pK*), but by 63 + 264 instead for ac(p). This difference has 
been discussed in [51]. In order to account for the corresponding S'f/(3) breaking, a 
rather generous correction factor k (168) has been allowed for. While this is certainly a 
valid procedure, it is somewhat against the spirit of using experimental data to constrain 
the penguin in B^, — > PlPI- An unexpectedly large penguin annihilation effect in this 
channel, beyond the available model estimates, would not necessarily be indicated by 
the B~ — s> K^pJ^ rate, not even in the SU{3) limit. In this respect, the method of 
[51] amounts to the standard analysis of CP violation in Ba PlPl with input on 
from factorization calculations, which are validated by comparing similar theory results 
on the penguin mode B~ — > K^p'l with data. Indeed, QCD factorization works very 
well for K*^Pl with little room for sizable power corrections. Correspondingly, 

the values for Tp and the angle 7 determined in Table 9 are very close to the values 
found in the factorization analysis, eqs. (135) and (143). Nevertheless, an independent 
control of penguin annihilation corrections in B^ — > PlPl^ which is not guaranteed by 
S~ — > K'^pI, would be very desirable. A variant of the method in [51] that can provide 
this control will be discussed in the following section. 
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5.6 Extracting Vp from Bd 



In this section we propose a method to constrain the penguin parameter rp in Bd — > Pi,Pl 
(133) using SU (3) flavour symmetry and data on the penguin decay Bd —>■ K^K}^. This 
approach shares the basic idea with the method discussed in section 5.5. An important 
difference is that now, unlike the case of section 5.5, the penguin process exhibits an 
exact SU{3) relation to the penguin amplitude of B^ PlPl- Because this relation 
extends beyond the heavy-quark limit, the method offers an independent control of all 
power corrections, in particular those from weak annihilation topologies. We show that 
a precise determination of the unitarity triangle is possible, already with present data 
on Bd KfKf. Since the penguin decay Bd KfKf is a A,S = transition, the 
up-quark sector of the amplitude does not have the same CKM suppression as for the 
AS* = 1 process B~ K*^ p~^. We will find that it is still sufficiently well constrained. 

The 5't/(3) relation between the relevant penguin amplitudes can be demonstrated 
as follows. The penguin contribution for Bd PlPl given by the component of 
the amphtude proportional to Ac = VcbV^- The corresponding part of the effective 
Hamiltonian (1) has the form 



where we have neglected higher-order clectrowcak effects. The operators Qi are defined 
in (2). The Hamiltonian in (171) gives rise to the QCD penguin amplitude in the charm 
sector of both Bd — > PlP2 and Bd — > K^K^. To prove the symmetry relation we 
note that all operators entering (171) are invariant under SU{2) rotations of the doublet 
(u.s) of quark flavours, the V-spin subgroup of flavour SU{3). The initial state, a Bd 
meson in both cases, is likewise a V-spin singlet. The final states p^p^ and K*^K*'^ 
are transformed into each other by interchanging u and s quarks, which represents a 
particular V-spin rotation. In the V-spin symmetry limit, therefore, the relation 



holds as an identity in QCD. As a consequence, the QCD penguin amplitudes propor- 
tional to Ac in Bd — * PlPl K'^K^ have the same form, including the weak 
annihilation contributions. This can be seen from (58), (59) and (95), (96). 

In practice V-spin is broken because the masses of up and strange quarks are not 
the same. This source of V-spin breaking can be expected to be of the typical size of 
flavour SU{3) breaking effects, roughly 20-30%. It is possible to estimate the required 
correction to the V-spin limit using factorization. We will give a more quantitative 
treatment below. Electroweak effects also violate V-spin symmetry. They are similar to 
isospin breaking and hkely to be much smaller than the 5'C/(3)-breaking effects due to 
the strange-quark mass. For example, the relative importance of (standard) electroweak 
penguins is governed by the ratio 010/^4 ~ 0.03. This is safely negligible in comparison 
with the dominant V-spin breaking effects. Contributions from electroweak penguin 
annihilation are very small and can also be neglected. 



"^QCDP.c = — ^ iCiQ'{ + C2 Q2 + ^ Ci Qi + Csg Qsg ] + h.c. 




(171) 



{Kl'Kl'\nQCBP,c\Bd) = {ptPl\nQCBF,cW 



(172) 
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We next turn to the phenomenological implications of the flavour symmetry relation 
(172). Denoting by ap{K*) = ap{Bd K*^K*^) the coefficient of {iGF/V^)\ in the 
amphtude for B^. — > K*^K*^, eqs. (59) and (96), the CP averaged branching fraction 
may be written as 



I \ 12 

*0 7>'*0\ _ ^Ba^F I'^c] 



327rmB 



foK{K*)\' + 2/i RealiK*) A{K*) + f2\A{K*)\' 



(173) 

A{K*) = a,{K*) - a^{K*) (174) 

where the functions fi depend only on CKM parameters. Expressed in terms of r = cot /3 
from (130) and 

a = cot 7 (175) 

they read 

-(^-^ 

In the region of interest for a Standard Model test, 

(7 = 0.447 ± 0.253, t = 2.54 ±0.13 (179) 

there is a clear hierarchy among the CKM factors 

/o = 0.835l°:ii, /i//o = 0.018 ± 0.086, h/ fo = 0.1611^^2 (^gQ) 

implying -C /o and /2 -C /o- The second inequality is a consequence of the fact that 
numerically |T4b/VtdP <^ 1. The first inequality arises because /i ~ cosa and the angle 
a is close to 90°. A similar feature holds for the decay B — > p7, where it leads to a 
suppression of hadronic uncertainties [52,53]. The dominance of the /o term in (173) is 
re-inforced by the hadronic factors since the difference |A(i^r*)| is systematically smaller 
than \ac{K*)\. This difference is a next-to-leading order effect in QCD factorization, 
whereas Oc is present at leading order. In addition, several terms cancel in the difference 
ttc — du- First, the NLO hard spectator corrections are identical in the c- and u-sector 
and drop out, which eliminates the uncertainty due to A^- Spectator effects in Oc — a„ 
can only come from penguin diagrams at NNLO ((9(q;^)), which are very small [54]. 
Second, also weak annihilation effects cancel in general, in particular those that are 
taken into account in our model estimate (coefficients hi). The only exception would 
be more complicated hd — > sdds annihilation topologies involving charm and up-quark 
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loops. These are both power and Zweig rule suppressed and not expected to give a 
significant contribution. Numerically we find 



\a,{K*)-a^{K*)\/GeY'^Omit',^^l{Ao) ^'oZ{a^) iSZ K) (181) 
Normahzed to the central value of |ac(-fr*)| = 0.084 we then have 

'"''^Hl,- = 0.25lg:!g (182) 

Together with the CKM factors from (180) we estimate a relative suppression of the third 
term in (173) by 0.010 (+0.012, —0.007) with respect to the first term. For the second 
term we estimate a relative size of at most 0.009 ± 0.043, neglecting the phase between 
Qc and A. In this case the potential magnitude of the correction depends strongly on 
the CKM suppression due to /i, which can be checked after the CKM factors have been 
determined at the end of the analysis. 

Because of the smallness of the /i, /2 terms, the first term in (173) determines the 
branching fraction to very good approximation. In the SU (3) limit, and up to negligible 
corrections from electroweak penguins, ac{K*) is equal to the penguin amplitude adp) 
in S — > pjp^ (in a corresponding normalization). Introducing the SU{3) factor 



adp) 

we obtain the ratio of CP averaged branching fractions 



1.28 (183) 



B{B ^ PIPl) p^ + r]^ + rl + 2prpC0S(j)p 



(184) 



This result constrains the penguin parameter in — > p^Pi- The four variables p, 
fj, Tp and (j)p may now be determined from the four measurements of r = {1 — p)/fi, 
Sp{p,fj,rp,(l>p), Cp{p,fj,rp,(l>p) and b = B{Ba ^ K*^Kf)/B{B ^ p+p£). This analysis 
then depends on a single theoretical parameter, the S'f/(3) factor |^| = 1.28 ± 0.14, 
where we adopt the estimate in (183) and assign a 50% error on the magnitude of S'f/(3) 
breaking. The quantity ^ = ac{K*) / ac{p) is real to very good approximation, ^ pa |^|. 

The expressions for S and C in terms of p, fj, r and (f) are identical to the case of 
Bd — > n'^TT' discussed in [19,20]. They read 

^ _ 2f][p^ + ff -r'^ - p{l - r^) + (p^ + ff - l)r cos 0] 
~ {{I- p)"^ + ff){p^ + ff + r"^ + 2rp cos 4>) 



g= ..^.." (186) 



2rf] sin - 
p^ + ff + r"^ + 2rp cos ( 

We remark that discrete ambiguities in the determination of p, fj, r and (f) can be avoided 
using other constraints on the unitarity triangle, which exclude p, fj far outside the region 
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allowed in the Standard Model [19,20]. A discrete ambiguity in the sign of cos0p can 
be resolved by the heavy-quark limit, which favours the solution with cos0p > 0. As 
pointed out in a similar context in [51], the discrete choice is still less restrictive in 
practice, because the second solution has cos^^ < —0.8, which is in fact far smaller than 
zero. 

The result of this analysis is given in Table 10, where we have also summarized the 
experimental input from Table 1 and section 5.4. The output values p, 77, Vp and (j)p are 

Table 10: CKM and penguin parameters extracted from r = cot /3 = 2.54 ± 
0.13, Sp = -0.05±0.17, Cp = -0.06±0.13 and b = B{Bd KfKf)/B{B 
PlPl) = 0.043 ± 0.015. The SU{3) breaking parameter is taken to be |^| = 
1.28 ±0.14. 





central 


T 


Sp 




b 




p 


0.088 


-0.010 
+0.011 


-0.029 
+0.029 


-0.000 
+0.010 


-0.009 
+0.011 


+0.006 
-0.007 


V 


0.359 


-0.014 
+0.015 


+0.011 
-0.011 


+0.000 

-0.004 


+0.004 
-0.004 


-0.002 
+0.003 


7[deg] 


76.2 


+1.0 
-1.0 


+4.7 
-4.9 


+0.0 
-1.6 


+1.5 
-1.8 


-0.9 
+1.1 


Tp 


0.064 


-0.003 
+0.004 


-0.002 
+0.003 


+0.000 
-0.000 


+0.010 
-0.012 


-0.006 
+0.008 




-0.20 


+0.00 
-0.00 


+0.00 
-0.00 


+0.43 
-0.47 


+0.03 
-0.05 


-0.02 
+0.02 



shown together with their sensitivity to the relevant input quantities. From Table 10 we 
draw the following conclusions: 

a) The errors on the CKM quantities p, fj and 7 are rather small. They are dominated 

by the uncertainty in Sp {fj is sensitive also to r = cot (3). 

b) The error from the SU{3) factor |^| is smaller than the errors from the experimental 

quantities r, Sp, b, which may still be improved by future measurements. 

c) The penguin parameter is obtained as = 0.064 ±0.014. The central value is some- 

what larger than the theoretical number in (135), but both results are compatible 
within errors. This confirms the expected smallness of r^, which is the basis for a 
precise extraction of CKM quantities. 

d) The phase (pp is seen to be strongly dependent on Cp, but essentially uncorrelated 

with the remaining parameters and input quantities. In particular p, f] and 7 are 
almost unaffected by the value of Cp within the measured range. This behaviour 
is in agreement with the general expectation discussed in 5.4. The error on (pp is 
completely dominated by the error on Cp. The sign of (f)p is opposite to the cen- 
tral standard model value from factorization in (136). If higher-order perturbative 
corrections cannot account for this change in sign, and we assume it is not due to 
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New Physics, this would mean that power corrections give an important contribu- 
tion to the strong phase. A similar situation is known to occur for the direct CP 
asymmetries in Ba — > tt'^tt' and Ba — > 7r'^K~. However, within uncertainties the 
numbers for (f)p in (136) and Table 10 are fully consistent with each other. The 
result for 0p in Table 10 confirms the prediction of a suppressed phase in the heavy 
quark limit. 

Combining the errors in Table 10 we find for the CKM angles 

7 = (76.2 ±5.3)° q; = TT - /3 - 7 = (82.3 ± 5.4)° (187) 

where the uncertainty is dominated by the experimental error in Sp. The result is in 
very good agreement with (143). It is already rather accurate at present. Prom Table 
10 we see that a precision of ±1° for 7 and a from this method should be possible. 

Pinally, we remark that the approximations leading to (184) may be cross-checked 
using the extracted value of 7 or o" = cot 7 = 0.25 ± 0.10, and r = cot j3 = 2.54 ± 0.13. 
Varying also the hadronic input parameters, the relative importance of the correction 
terms in (173) is then smaller than ±3%. (Por the default parameter set and a — 0.25, 
T = 2.54, the correction is —0.8%.) The corresponding change in (184) could be absorbed 
in a modification of |^| by less than ±1.5%, which is entirely neghgible. 

5.7 Unitarity triangle from Bd — > and Bd — > K^K^ 

The analysis of section 5.6 made use of CP violation in B^ ^ PlPl^ ^ nieasurement of 
sin 2/9, and the B^ K*^K*^ branching fraction to obtain an accurate determination of 
the unitarity triangle. The decay B^ — > K*^K*^ served to fix the penguin-to-tree ratio 
in Bd —>■ PlPl based on SU{3) flavour symmetry. 

The same analysis may also be performed with the y^-modes replaced by their 
pseudoscalar counterparts, that is, using CP violation in B^ — > tt'^tt^ and constraining 
the penguin parameter r^^ with B^ K^K^ and SU{?>) symmetry. The formulas of 
section 5.6 apply with obvious substitutions. Related discussions can be found in [19,20]. 

Using form factor estimates based on [55] 

/:^^'^(0) = 0.258 ± 0.031 /+^-^(0) = 0.304 ± 0.042 (188) 

we find from QCD factorization [6,23] 

= 1.46 ±0.23 (189) 

Again we have assigned a generous 50% uncertainty on the total amount of S'f/(3) 
breaking. With experimental input from [24] we obtain the results displayed in Table 
11. 

Combining errors we obtain from Table 11 

7 = (61.1 ±6.8)° (190) 



ICI 



ac{K) 
ac(vr) 
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Table 11: CKM and penguin parameters extracted from r = cot /5 = 2.54 ± 
0.13, = -0.61 ± 0.08, = -0.38 ± 0.07 and b = B{Bd K^K^)/B{B 
7r"'"7r~) = 0.186 ± 0.040. The SUii) breaking parameter is taken to be \C,\ — 
1.46 ±0.23. 





central 


r 






b 


ICI 


p 


0.179 


-0.013 
+0.014 


-0.020 
+0.023 


-0.014 
+0.021 


-0.015 
+0.019 


+0.023 
-0.027 


V 


0.323 


-0.011 
+0.011 


+0.008 
-0.009 


+0.005 
-0.008 


+0.006 
-0.007 


-0.009 
+0.011 


7[deg] 


61.1 


+1.0 
-1.0 


+3.3 
-3.7 


+2.3 
-3.4 


+2.5 
-3.1 


-3.7 
+4.4 




0.146 


-0.008 
+0.009 


-0.006 
+0.008 


-0.002 
+0.002 


+0.015 
-0.017 


-0.020 
+0.027 


07r 


-0.88 


+0.01 
-0.01 


+0.02 
-0.03 


+0.19 
-0.24 


+0.09 
-0.13 


-0.17 
+0.14 



This value is lower than the result in (187) but it remains consistent at the level of roughly 
2a. One possible source of this discrepancy is the rather large value of C^r = — 0.38±0.07, 
representing the simple average of the BaBar [56] and Belle [57] results 

= -0.21 ± 0.09 (BaBar) = -0.55 ± 0.09 (Belle) (191) 

These results are not in very good agreement. With a smaller ICj^l, favoured by QCD 
factorization and the BaBar measurement, the extracted value for 7 would increase 
somewhat. For example, with Ctt = —0.1 we obtain 7 = 66.6°, keeping all other inputs 
fixed. Particularly important for the resulting 7 is the value of St^. If it were 2a lower 
in absolute magnitude, at 3.,^ — —0.45, the central value of 7 would shift to 7 = 67.4°. 
The uncertainties in b and \(\ also have a relatively large impact. This is because of the 
larger size of the penguin contribution r^^ in comparison with Vp. Note that the error on 
7 from the uncertainty in b is almost twice as large in Table 11 than in Table 10, even 
though b is known with an accuracy of 22% in the former case and only to 35% in the 
latter. Thus, because of the larger size of the penguin amphtude, and also because of 
the experimental situation of C^, which is still not entirely resolved, the determination 
of the unitarity triangle from B tt+tt" and B K^K^ appears to be somewhat less 
precise than the determination from B — > p^p^ and B — > K'^K'^. 

5.8 CP violation in Bs — > <Pl<Pl 

The decay Bg is a pure b ^ s penguin transition and thus of considerable interest 

as a New Physics probe. Possible hints of deviations from the Standard Model in CP 
violation in the 6 — > s penguin process B^ —>■ (j)Ks, and similar modes, have so far 
remained inconclusive. A detailed experimental study of Bg — > 00 will become possible 
with the LHC [1]. In the Standard Model CP violation in Bs (jxp is small. Any 
nonzero effect in excess of the Standard Model contribution will signal the presence of 
New Physics. Based on our next-to-leading order results we shall investigate the size 
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and uncertainty of CP violation in the Standard Model, which ultimately limits the 
sensitivity to New Physics. 

The time dependent CP asymmetry in Bg ^ (f)L<^L decays is defined by 

T{Bs{t) (Pl<Pl) + ^{Bs{t) (Pl<Pl) 

(192) 

Here we have neglected the effects of a nonzero width difference AF^^, which would 
modify the time dependence of the CP asymmetry. This can be taken into account in 
extracting and C^, but would not change the following discussion of these parameters. 
For a generic B decay into a CP self-conjugate final state / one has 

c_ 21me ^_ , _ Ml, A{B /) 

where = {B\Hab=2\B) is the B-B mixing amplitude. We use the phase convention 
CP\B) — — \B). The CP violation parameters -S",^ and then become 

S^^2X T] Re T-T , C^^2X T] Im —r (194) 

ac(0) ac((/)) 

where ap((f)), p = u, c, is the coefficient of {iGF/^/2)Xp in the Bs (j)L4'L amplitude (74) 
and (113). It can be seen from (194) that and depend on the same CKM quantity 
but on different hadronic parameters. A measurement of is therefore only of limited 
use in controlling hadronic uncertainties in 5"^. 

The hadronic parameters in (194) depend on the difference between the penguin 
amplitudes from the charm and the up-quark sector. This difference is calculable in 
factorization. It has the further advantage that the leading annihilation corrections 
related to the parameters bi (113) cancel in ac{(p) — au{(p)- A similar cancellation occurs 
for the hard-spectator scattering contributions in the NLO approximation. We then find 

|ac(0)-a.(0)|/GeV3 = O.O57lHS^(^) («2^) t'ofiU^c) ;2;o2i (/^) (195) 

where we show the dominant parametric uncertainties and their origin (in brackets). 
In contrast to the difference ac{(j)) — aii(0), the absolute value of ac(0) depends on the 
annihilation contributions 6j. Rather than aiming for an accurate theoretical prediction, 
it therefore appears more reliable to extract |ac(0)| from experiment. Neglecting the 
very small up-quark contribution, we may write 

647rmB^ 

This gives 

(197) 



0.177 GeV^ 


"S(5, -^0i0i)" 


1/2 


1.53 ps 


1/2 




15 • 10-6 
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First evidence for the decay Bg — > (jxp has been reported by the CDF collaboration, 
which quotes [58] 

B{Bs #) = (Uttistat) ± 6{syst)) ■ 10"^ (198) 

This is not yet the longitudinal branching fraction needed in (197), and the error is 
still large. By the time CP violation in Bs — > (fxp will be studied at the LHC, the 
branching fraction will be known with good precision and the number in (197) can be 
easily updated. 

The quantity is predicted to be small and positive in the Standard Model. With 
our default parameter set and 77 = 0.36 we obtain 

5<^(default) ^ 0.01 (199) 



From the discussion above we conclude that the Standard Model upper limit can be 
written as 

^ 2X^r, ^ X^r, (200) 

\ac{(p)\ 



15 • 10-6 

A similar limit holds for the absolute value of C,^. For rj < 0.4 we have 

S4, ^ 0.02 \C^\ ^ 0.02 (201) 

A rescaling for the actual value of the branching fraction, which might deviate from the 
assumed default value of 15 • 10"^, can be done using (200). Measurements in excess of 
these Standard Model limits would constitute evidence for New Physics. The expected 
sensitivity of LHCb after five years of data taking is cr(5'0) ~ 0.05 [59]. Improvements 
to values of 0.01 or 0.02 with the anticipated LHCb upgrade appear possible [1]. This 
should allow us to exploit the full New Physics potential of and to detect non-standard 
effects in 6 — > s penguins at the few percent level. 



6 Comparison with the Hterature 

In this section we comment briefly on the existing literature related to the subject of the 
present paper [9-18]. 

We re-emphasize that factorization calculations for charmless two-body B decays, in 
particular B — > VlVl, are useful for flavour physics analyses such as the determination 
of CKM parameters. This has also been stressed in [11] and earlier in [6,19,20,23]. 
The most comprehensive study of B ^ VV decays has been presented in [11] with 
an emphasis on total branching fractions and polarization observables, for instance the 
longitudinal polarization fractions /l. More recently these processes were considered in 
[18] in an extended study, following the analysis of [11]. We do not discuss transverse 
polarization here but rather concentrate on the decays with longitudinal vector mesons 
B VlVl. These arc calculable in QCD in the heavy-quark limit and thus of special 
interest for phenomenological applications in flavour physics. We list the detailed results 
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for the B — > VlVl amplitudes in explicit terms. The corresponding results of [11] can 
be reconstructed from similar formulas given for B — > PV decays in [23]. Our main 
results are consistent with [11]. A minor difference with (the original version of) [11] are 
the expressions for penguin annihilation A\ ^ and (84). The final expressions in 
[11] give incorrectly, due to a relative sign change, a nonvanishing and A^ ^ 0, even 
though the basic formulas agree with (82). The difference leads to a reduced sensitivity 
to penguin annihilation in penguin-dominated B — > VlVl decays. It does not play a role 
for B — > p" 1 because the corresponding decay amplitude is tree dominated and color 
allowed, so no deviations from modelling of power-suppressed contributions are expected. 
This point has previously been noted in [18]. There is consensus on the expressions (82), 
(84) within the annihilation model used [9,18,61] . Another difference with [11,23] is 
the treatment of electromagnetic penguin matrix elements contributing to a^ g, where 
we have proposed an explicit model for the long-distance contributions to these 0{a) 
terms. 

The article [9] concentrates on transverse polarization and therefore docs not report 
the complete expressions for the amplitudes with longitudinal vector mesons. Where a 
comparison is possible we agree with the results of [9], except for two minor discrepancies. 
One is the detailed form of the integrand of the annihilation parameter A^^. However, 
the final result for A{^ coincides with ours. Another difference is the annihilation part 
of B~ — > K*^(t), which should read 63 -|- instead of 63 — 6f^/2. Both issues are 
inconsequential. 

B — > VV decays have been studied within QCD factorization also in [13-18]. These 
papers address various VV channels, especially penguin dominated modes such as ^X*. 
Some of them investigate the impact of New Physics scenarios [15,16,17], [18] extends 
the analysis to V A and AA modes as well. In comparison, the present paper, while 
concentrating on B ^ VlVl, gives complete NLO results for all channels, a detailed 
analysis of uncertainties and applications for precision tests of flavour physics. The 
authors of [18] employ mc{mb) — 0.91 GeV, smaller than the value used here and in [11]. 
We find that the error due to rric for the longitudinal amplitude is small compared to 
other experimental input, also for penguin dominated decays. 

7 Conclusions 

In this paper we presented a systematic analysis of i?-meson decays into a pair of longi- 
tudinal vector mesons. The main results can be summarized as follows: 

• Explicit formulas are given for the complete set of AS = and AS = 1 decay 
amplitudes of 5 — > VlVl at NLO in QCD factorization. Estimates of power 
corrections from weak annihilation are included to study the sensitivity to effects of 
this kind in phenomenological applications. The set of decays considered comprises 
17 AS* = and 17 AS — 1 channels, including 2 and 4 pure annihilation modes, 
respectively. 
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The agreement with the available measured branching ratios of Bd, B decays 
into p+p-, pV°, P"P°, K*^K*^, K*^P^, K*-p+, K*-p^, K*^P', K*y, K*-(j) and 
K*^uj is very good, within current uncertainties and with the central values used 
for B ^ Vl form factors. The decay to p~uj agrees at the level of about 2a. We 
note that QCD factorization works well in particular for the penguin modes and 
for the three pp channels and their characteristic hierarchy of branching fractions. 
Our hadronic input is based on the available literature. No tuning of parameters 
has been done to improve the fit with data. 

Long-distance electromagnetic penguins in B decays with neutral vector mesons 
p^, cu or (p are taken into account using a model description. They are of some 
conceptual interest but their numerical impact is generally small. 

The deviation from ideal mixing in the a;-0 system is found to have a small effect 
on most decay modes with these particles in the final state. The impact is very 
large for B~ — > p~(f). 

We advocate the direct use of mixing- induced CP violation in B^ — > P^Pl, mea- 
sured by Sp, to extract the parameters of the unitarity triangle. Together with 
sin 2/3 the current measurements of Sp imply (143) 

7= (72.4 ±6.2)° (202) 

where the error is dominated by Sp. This analysis benefits from the small penguin- 
to-tree ratio for vector modes Vp — 0.038 ± 0.024 (135), which leads to a residual 
theory uncertainty in 7 of ±3°. 

We propose a method to relate the penguin contribution in Bd — > P^Pl ^'^ 
decay Bd K*j^Kf based on the V -spin subgroup of flavour SU{'i). This makes it 
possible to constrain the uncertainties due to penguin power corrections, especially 
from annihilation topologies, and provides us with a check on the penguin-to-tree 
ratio calculated in QCD fatorization. The absolute value determined by the F-spin 
method, Vp — 0.064±0.014, is consistent with the calculation in QCD factorization. 
The resulting angle 

7 = (76.2 ±5.3)° (203) 
has a residual theory error of ±1°. 

A comparison of the analyses mentioned in the previous two items to the corre- 
sponding ones with pseudoscalar decay modes suggests that the hadronic uncer- 
tainties are better under control in the case of vector modes. 

We point out that within the SM sin 2/3 and the CP violation parameter Sp in 
B ^ P^P2 determine 

iV^bl = (3.54 ± 0.17) • 10~^ (204) 
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where the error is at present still entirely dominated by sin 2/3. Hadronic uncertain- 
ties enter only at second order in Sp and the penguin parameter and are below 
2%. Possible New Physics affecting the Bd — mixing phase can be constrained 
by comparing the above value of \Vuh\ with direct determinations from exclusive or 
inclusive h ^ ulf decays. 

• In future measurements Bs (jxf) will provide tests for New Physics. We present 
a bound on S^p and C^, which will further improve when more data are available. 

The phenomenology of i? ^ VlVl decays is rich and promising. QCD factorization 
provides a solid theoretical basis for these processes. With the measurement of additional 
channels and improved precision for the ones already observed, many more applications 
in flavour physics may be foreseen. 



A Long-distance electromagnetic penguins 

For B decays with neutral vector mesons p°, uj or (p, the operators Q\ 2 have electromag- 
netic penguin-type matrix elements where the photon from the pp loop is transformed 
into one of these mesons. The photon virtuality = my is then small. While the pen- 
guin loop may still be considered short-distance dominated for p = c due to the charm- 
quark mass, the matrix element becomes sensitive to long-distance hadronic physics for 
p = u. This can be seen from the perturbative result for the case of the charm quark 
(38), which diverges in the limit — > 0. The up-quark contribution is thus not strictly 
calculable. Since this situation arises only in a small electromagnetic correction, it is 
not a serious problem for practical purposes. In fact, additional dependence on long- 
distance hadronic physics is to be expected when electromagnetic radiative corrections 
to hadronic B decays are considered. Still the penguin matrix element under discussion 
contributes within our approximation scheme of including leading electroweak effects. 
We shall therefore give an estimate of its size using available information on the long- 
distance dynamics of the up-quark loop. Apart from obtaining a numerical evaluation 
of the effect, the long-distance electromagnetic penguin is also interesting for conceptual 
reasons. 

The up-quark loop is closely related to the vacuum polarization function Il{k'^), where 
the UV subtraction is given by the standard renormalization prescription of the weak 
hamiltonian. We thus write the penguin matrix element, needed at low photon virtuality 

— my, as the matrix element evaluated at k^ — ml plus a remainder proportional to 
the difference n(A;^) — Il{mf). The contribution of the electromagnetic up-quark penguin 
to the coefficients g then takes the form 

The correlator 11 (/c^) is defined through 




(0|Tj^(a;)j,(0)|0) = {k^k, - k^g.^Mk^) (206) 



{Ii{mi)-Ii{mi)) 



(205) 
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where = ifJuU. The first terms in the square brackets of (205) come from the per- 
turbative evaluation of the matrix element at k'^ — ml and carry the appropriate scale 
and scheme dependence. The remainder depending on U may be computed to lowest 
(one-loop) order, which reproduces the perturbative result for the matrix element. We 
shall treat n(A;^) — Il{ml) as the full hadronic correlator, which includes the nonper- 
turbative hadronic physics relevant at low k"^. This procedure assumes a factorization 
of the soft hadronic correlator form the remaining parts of the diagram, which is not 
strictly justified. We adopt this additional assumption to obtain a rough estimate of the 
long-distance sensitive penguin contributon. A similar method has been proposed and 
applied in the context of b ^ s{d)e^e~ decays in [60]. 
The function Il{k^) obeys the dispersion relation 

ny) = ir,i« . (207) 

^ ^ n Jo t-k^ -le ^ ' 

In this form the dispersion relation needs one subtraction, but the subtraction constant 
cancels in n(A;^) — n(m^). In principle Imn(t) could be determined experimentally. 
Instead, for simplicity, we choose a convenient ansatz that should capture the essential 
features of the true hadronic quantity Imn(f). We write Imll as the sum of a resonance 
and a continuum contribution 



Imn = Imn^ + ImHe (208) 



where 



r=p,ijj ^ ' ' ' ' 

Imne(t) = ^0(tc -t) + ^e(t - Q, = ATT^f^ +fl), t>0 (210) 

The asymptotic QCD result fixes Imllc to Nc/ (127r) = 1/ (An) at large t. Imposing quark- 
hadron duality for the integral of Imn(t) up to (at least) t = determines the value of 
tc = Air'^^fp + /J) ?a 3.1 GeV^. The factor 1/2 in (209) is an isospin factor coming from 
the overlap of p° and cu with the ■u7^m current. Determining 11 in (205) with the help 
of (207) and (208), treating the resonances as narrow and taking the heavy-quark limit 
tc ^ rrLf,, we finally obtain (39). Concerning the factor in square brackets in (39), two 
limiting cases are worth noting. If k"^ = rriy — ^ 0, we recover an expression similar to 
(38) where the light-quark mass under the logarithm is replaced by the hadronic scale 
^/tl. In the limit k'^ — rriy — > tc the same terms appear, and in addition the perturbative 
imaginary part —2m/?,. 



B Coefficients a^, hi 

In the following Table we quote the central values of the coefficients at as defined in 
(23) for two final-state p-mesons (3). The default value used for the model of power- 
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suppressed hard-spectator contributions is Xh = In ^ and for the renormahzation scale 
it is // = 4.2GeV. 



ai 


02 


03 + 05 




04 


0.991 + 0.020/ 


0.177 - 0.081/ 


0.002 - 0.001/ 


-0.025 - 0.016/ 


-0.033 - 0.009i 


(«.i' + «;;)/n 








"10/' » 


-1.84 - 0.54i 


1.15 + 0.02i 


-1.10-0.02i 


-0.17 + 0.09i 


-0.17 + 0.09i 



The central values of the coefficients hi as defined in (78) for two final-state p-mesons 
(3) are given below. The default value used for the model of power-suppressed annihi- 
lation contributions is Xa = In ^ and for the renormahzation scale it is ji — 4.2GeV. 
Here ta = Bp^jApp. 



TAbi 


rA&2 




rA&4 


TAbi'^/a 


TAbf'^/a 


0.029 


-0.011 


0.003 


-0.003 


-0.035 


0.013 
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